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FOREWORD 


During the early part of World War II, some of the helicopters 
designed for military use were observed during ground tests to exhibit 
a violent oscillatory rotor instability which endangered the safety of 
the aircraft. This instability was at first attributed to rotor -blade 
flutter, but a careful analysis indicated it to be caused by a hitherto 
unknown phenomenon in which the rotational energy of the rotor was con- 
verted into oscillatory energy of the blades. This phenomenon was 
usually critical when the helicopter was operating on or near the ground 
and, hence, was called ground resonance. An oscillatory instability of 
such magnitude as resulted from this phenomenon would generate forces 
that could quickly destroy a helicopter. The research efforts of the 
National Advisory Committee for Aeronautics were therefore enlisted to 
investigate the difficulties introduced by this phenomenon. During the 
interval between 1942 and 1947, a theory of the self -excited instability 
of hinged rotor blades was worked out by Robert P. Coleman and Arnold M. 
Feingold at the Langley Aeronautical Laboratory. This theory defined 
the important parameters and provided design information which made it 
possible to eliminate this type of instability. These results were 
originally released by the NACA in three separate papers, as follows: 

(1) NACA Advance Restricted Report 3G29, 1943 (Wartime 

Report L-308) entitled "Theory of Self -Excited 
Mechanical Oscillations of Hinged Rotor Blades," 
by Robert P. Coleman 

(2) NACA Advance Restricted Report 3H3> 1943 (Wartime 

Report L-312) entitled "Theory of Mechanical 
Oscillations of Rotors With Two Hinged Blades," 
by Arnold M. Feingold 

(3) NACA Technical Note 1184, 1947, entitled "Theory of 

Ground Vibrations of a Two-Blade Helicopter Rotor 
on Anisotropic Flexible Supports," by Robert P. Coleman 
and Arnold M. Feingold 

These three reports present the fundamental reference material on the 
subject of rotor mechanical instability but have been for some time out 
of print. As a result of demands for this type of information, they are 
being combined in the present volume and reissued with appropriate correc- 
tions. The combined volume consists of three chapters representing the 
three aforementioned papers in consecutive order. 

Inasmuch as the authors of these papers are no longer with the NACA, 
the task of checking these papers and incorporating such corrections as 
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seemed proper was undertaken by George W. Brooks, Vibration and Flutter 
Branch, Dynamic Loads Division, of the Langley Aeronautical Laboratory. 
Mr. Brooks also prepared an appendix to chapter XXL, which deals with 
the general equations of motion for a two-blade rotor and includes the 
effects of damping. 


Arthur A. Regier 

Head, Vibration and Flutter Branch, 
Dynamic Loads Division, 
Langley Aeronautical Laboratory 
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THEORY OF SELF-EXCITED MECHANICAL OSCILLATIONS 
OF HELICOPTER ROTORS WITH HINGED BLADES 
CHAPTER I 

THEORY OF SELF -EXCITED MECHANICAL OSCILLATIONS 
OF HINGED ROTOR BLADES 
By Robert P. Coleman 

SUMMARY 


Vibrations of rotary-wing aircraft may derive their energy from the 
rotation of the rotor rather than from the air forces. A theoretical 
analysis of these vibrations is described and methods for its application 
are explained herein. 

The theory includes the effects of unequal stiffness of the pylon 
for deflections in different directions and the effect of damping in the 
hinges and in the pylon. Both the derivation of the characteristics equa 
tion and the methods of application of the theory are given. In particu- 
lar, the theory predicts the so-called "odd frequency" self -excited speed 
range as well as the shaft critical speed. Charts are presented from • 
which the shaft critical and the self -excited instabilities can be pre- 
dicted for a great variety of cases. The influence of each physical 
parameter upon the instabilities has been obtained. The comprehensive 
treatment applies to a rotor that has any number of blades greater than 
two. O nl y a brief discussion and the formula for shaft critical speed 
are given for the one- or two -blade rotor. 

The use of complex variables in conjunction with Lagrange's equa- 
tions has been found very convenient for the treatment of vibrations of 
rotating systems . 


INTRODUCTION 


A rotary -wing aircraft that has hinged blades will, under certain 
conditions, be subject to vibrations which derive their energy from the 
rotation of the rotor instead of from the air forces . The term "ground 
resonance" usually refers to vibrations of this type. Although such 
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vibrations have apparently caused accidents in some rotary-wing aircraft 
and have impaired the flying qualities of others, very little attention 
has been given this problem in the literature. A theoretical analysis 
has therefore been undertaken, and the purpose of the present paper is 
to present the theory and to describe the application of the theory to 
rotary -wing aircraft. 

General vibration theory and its application to allied problems as 
well as to the particular problem of rotor vibration are discussed in 
references 1 to 4. A good general background for the present problem is 
provided in the chapters on rotating machinery and on self -excited vibra- 
tions in reference 1. References 2 and 5 treat in more abstract fashion 
the topics of rotation and damping. A discussion of the variety of modes 
of vibration that exist in rotors and. a number of frequency formulas 
obtained by considering separately each degree of freedom are given in 
reference 4. This discussion does not, however, lead to a prediction of 
self -excited modes of vibration. 

Experience has shown that two types of mechanical vibration may 
occur in rotors. The vibration frequency of the pylon is equal to the 
rotational speed in one type and is unequal in the other. The first 
type is sometimes called the even-frequency vibration or the one-to-one 
frequency, and the second type, the odd frequency. The one-to-one fre- 
quency vibration resembles the phenomenon occurring at a critical speed 
of the shaft of rotating machinery and will consequently be referred to 
in this paper as a shaft critical vibration. The odd -frequency vibra- 
tion is properly called a self -excited vibration. 

A derivation of the characteristic equation for the whirling speeds 
of a three -blade rotor has been given by Wagner of the Kellett Autogiro 
Corporation. By considering only the case of a pylon having equal stiff- 
ness in all directions of deflection, Wagner has shortened the analysis 
by considering directly the equilibrium of forces and moments under con- 
ditions of steady circular whirling. Some examples of the dependence of 
whirling speed upon rotational speed are given, and the formula for the 
shaft critical speed is obtained. 

In the present paper, the theory also includes the effects of damping 
in the hinges and in the hub and the effects of different stiffnesses of 
pylon deflection in different directions. The method of analysis, partic- 
ularly the use of complex variables in the equations of motion, is explained 
in some detail and all the previous results are shown to be a special case 
of the more general problem treated here. 
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SYMBOLS 


radial position of vertical hinge 

All 
All 

Mu = ^11 1 

12 = A^2 I elements of the characteristic determinant (see eq. ( 31 )) 

a 21 = A 21 

A 22 . 

A 22 

h 
B 


distance from vertical hinge to center of mass of blade 
damping force per unit velocity of pylon displacement 

( Bx 

Bf = 


■P) 


ABf = 

B I 

% 

e, Cp, 

Cl 

CR 

•D 

F 


Bx “ By 


coefficient defined in equation (35) 
coefficient defined in equation (34) 
. arbitrary constants 

coefficient defined in equation (35) 
coefficient defined in equation (34) 

time-derivative operator, d/dt 
dissipation function 
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moment of inertia of blade about hinge, n^b^h. + 


I 1# ... I 5 coefficients defined in expressions ( 37 ) 

j, k indices and subscripts used with hinge coordinates (eq. (l4) ) 

(% - 


K 


spring constant 


AK f = 


Kx - % 


m x + 


m 


effective mass of pylon l mf = 


“b 


blade mass 


Amf 


mx - my 
2 


M total effective mass of blades and pylon, mg + nm^ 

AM mass added at hub for vibration test 

n total number of blades 

r radius of gyration of blade about its center of mass 

R]_, ... coefficients defined in expressions ( 37 ) 

s stiffness ratio, Ky/K x 


t 

T 

Tr 

Tk 

T s 


time 

kinetic energy 

kinetic energy of rotation of blade about its center of mass 
kinetic energy of translational motion of kth blade 
kinetic energy of pylon 


V 


potential energy 
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x, y displacements 

x 0 , yo values of x and y when t = 0 

z complex displacement, x + iy 

z conplex conjugate of z, x - iy 

a angle between blades, 2it/n 

Pq, Pi, ... P]j angular displacements of blades 

Pj^ value of 3^ when t = 0 

£q, ... variables representing hinge deflections when equations 

are expressed in fixed coordinate system 

9 0 , 0i, ... 0^ variables representing hinge deflections when equations 

are expressed in rotating coordinate system 

B x / Bv \ 

A-j, = — in applications 

X M x W ) 


\ — in applications ] 

V My \MyO>r / 

Ap = in applicationsj 


. Bf 
= M" 


■\ _ ®a 

Aa ' ~W 


AAf = 


ABp 


Al = 


b 1 + ^ 


A?= ^ 

in applications) 

I \lpr 2 ) 
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Ax = 


4 * S) 


mass ratio , 


nnib 


v±, v 2 

CD 

(Of 


nxp + nm^ 

expressions defined in equation ( 3 ) 

angular velocity of rotor (the dimensionless ratio co/a^ is 
called (jd in applications) 

angular whirling velocity measured in rotating coordinate 
system (used in nondimensional form in applications) 

angular whirling velocity measured in fixed coordinate system 
(used in nondimensional form in applications) 


< 1 ^ reference frequency, ^xMy. 

Subscripts : 


f fixed coordinate system 

a rotating coordinate system 

P hinge deflection 

x, y component directions in fixed coordinate system 


APPROACH TO THE VIBRATION PROBLEM 


Stability and Instability 


If a vibrator were attached to a rotorcraft, several modes of vibra- 
tion could be excited at any rotor speed. Only the modes that are likely 
to be excited during operation of the aircraft, however, are important. 

In the present discussion, it is convenient to classify the modes 
of vibration according to the circumstances required for their excita- 
tion. The different types of vibration are identified analytically by 
the nature of the roots of the characteristic equation. A hinged rotor 
may encounter three types of vibration which, for convenience, are herein 
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designated ordinary, self excited, and shaft critical. At zero or slov 
rotational speeds, an external force is required to excite vibration. 

The friction always present in such systems causes the vibration to be 
damped out when the force is removed. Modes of vibration requiring an 
external applied force to maintain them will be called ordinary vibra- 
tions . The mathematically idealized case of zero damping will also be 
considered an ordinary vibration when it is understood to be an approxi- 
mation to a system actually damped. Self -excited modes of vibration are 
those with negative damping and are recognized analytically by the fact 
that a root of the characteristic equation is a complex number which has 
a negative imaginary part. A slight disturbance will tend to increase 
with time instead of damping out. 

When a rotating system is not perfectly balanced, the centrifugal 
force of the unbalanced mass may excite vibrations that have peak ampli- 
tudes at certain rotational speeds . Vibration excited by imbalance and 
in resonance with the rotation will be called shaft critical vibration. 
This type occurs at the rotational speed at which the spring stiffness 
of the pylon is neutralized by the centrifugal force. In the analysis, 
the shaft critical vibration is recognized in rotating coordinates as a 
zero frequency and in fixed coordinates as a frequency equal to the 
rotational speed. The critical speeds of a rotating shaft are a common 
example of this class . 

The purpose of a theory of .rotor vibration is mainly to predict the 
occurrence of and, if possible, to show how to avoid self -excited and 
shaft critical vibrations . 


Choice of Degrees of Freedom 

Of the large number of degrees of freedom of a hinged rotor, the 
important ones for the present problem have been found to be hinge deflec- 
tion of the blades in the plane of rotation and horizontal deflections of 
the pylon. Other degrees of freedom, such as the flapping hinge motion 
of the blades, the bending or torsion of the blades, and the torsion of 
the drive shaft,- are considered unimportant in the problem of self -excited 
oscillations. Some motions, such as landing-gear deflection, that pro- 
duce lateral deflection at the top of the pylon may, however, be important. 


Physical Parameters 

The theoretical results given later provide a means of predicting 
the natural frequencies and, in particular, the critical speeds and 
unstable speed ranges in terms of certain physical parameters, such as 
mass, stiffness, and length. The successful application of the theory 
depends upon the determination of the proper values of these physical 
parameters for the aircraft or model being studied. 
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The important parameters to be determined are: 
a radial position of vertical hinge 

b distance from vertical hinge to center of mass of blade 

mass of blade (Flexibility of the blade structure may have to 
be taken into account by the use of an effective value of 
1115 different from the actual blade mass . The effective 
blade mass can be taken as the value required to make the 
theory predict the correct pylon natural frequency when the 
rotor has a zero or very slow rotational speed.) 

I moment of inertia of blade about hinge, mt>b 2 (l + 



Kp spring constant of self-centering springs, which can be deter- 

mined by a force test or from the hinge frequency with the 
hub rigidly supported 

m x , my effective mass of pylon for deflections in x- and y-directions 
K x , Ky effective stiffness of pylon 

The effective mass of the pylon is the value of a concentrated mass 
that would have the same kinetic energy expressed in terms of the deflec- 
tions at the hub as the actual pylon and hub if it were placed at the 
rotor hub in the plane of rotation. The effective stiffness of the pylon 
is the stiffness of a spring that, if placed at the hub in the plane of 
rotation, would have the same potential energy in terms of deflections 
at the hub as the actual pylon. Equivalent definitions are that, if a 
simple spring and mass were imagined to be substituted at the hub in the 
plane of rotation for the pylon and hub, . the natural frequency and the 
change of natural frequency with added mass would be the same as for the 
actual pylon . 

An experimental method of measuring the effective mass m x and 
stiffness K x of the pylon is to replace the rotor by an approximately 
equal, rigid, concentrated mass at the hub and to measure the nat- 

ural frequency for two or more values of this added mass. The quantities 
are then related by the equation 
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or 


-^2 = ^r - ( n bc + ^M) 

CDx . 

If measured values of l/ao* 2 are plotted against added mass £M and 
a straight line is drawn through the points, the intercept and the slope 
of the line will determine the effective values of K x and 1 % . 

For the parameters a and b, the actual geometric lengths should be 
used unless the flexibility of the hinge offset arm a is comparable in 
magnitude with the hinge spring stiffness. In this case, it is recom- 
mended that an effective value of a be guessed and that b be deter- 
mined by subtraction from the correct geometric value of a + b. 

The damping parameters may be defined by the form of a dissipation 
function F. The function 


2F = BxX f 2 + Byf/ + E p p k 2 

k=0 


is equal to the rate of dissipation of energy by damping. The param- 
eters Bx and By thus measure the damping force per unit velocity 

referred to linear displacements of the top of the pylon and Bp is 

the damping torque per unit angular velocity at a blade hinge . If the 
damping force per unit velocity is not a constant, effective values 
should be used that will represent the same dissipation of energy per 
cycle as actually occurs with a reasonable amplitude, of vibration. The 
amplitude of free vibration in a single degree of freedom is given in 
terms of Bx, By, and Bp by 


x f = x Q e 


ik t 

2M 


yf = y Q e 


2M 


21 


P k “ Pko e 
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The damping parameters are probably the most difficult ones to 
measure accurately . In practice, it is advisable to make calculations 
for a given case, first on the basis of no damping and then with the 
use of the estimated values of the damping parameters. 


MATHEMATICAL DEVELOPMENTS 


Method of Analysis 

The derivation of the characteristic equation that is used as the 
basis for predicting the unstable oscillations of a rotor is presented 
in this section of the paper. Readers interested solely in applications 
can omit this section and proceed immediately to the section entitled 
"Method of Applying Theory." 

The method of analysis treats the equations of motion for small 
displacements from the equilibrium condition with steady rotation. A 
proper choice of coordinates leads to equations with constant coeffi- 
cients . The solutions are exponential or trigonometric functions . 

The mathematical manipulations involved in treating the motions of 
a mass in a plane of rotation are facilitated by the use of a complex 
variable. The typical disturbed motion obtained by solving the equa- 
tions of motion is an elliptic whirling motion, which is represented in 
terms of a complex variable z = x + iy. At any instant, z represents 
the displacement of the pylon from its equilibrium position. An expres- 
sion such as 


iooft 

z = ce J - 

represents whirling of the pylon in a circle of radius c with angular 
velocity cop. The sign of o>f determines the sense of the rotation. 
Two rotations in opposite sense with the same radius are equivalent to 
a vibratory motion in a straight line and are given as follows : 

3 . c (e 1 ^ + e- 1 ^) 

= 2c cos (JOft 


Two opposite rotations of different radii- are equivalent to whirling in 
an ellipse. A complex value of represents whirling in a spiral, 

which may he either a damped or a self -excited motion depending upon 
the sign of the imaginary part. A self -excited motion exists when the 
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imaginary part of (Df is negative, and the magnitude of z increases 
with time. 

The displacements may he referred to a fixed or to a rotating coor- 
dinate system. If Zf and z a are the displacements with respect to 
a fixed and to a rotating reference system, respectively, 

Zf = Z a e lc0t 
z a = ce i«at 

then 

z f = ce i ^‘ KD ) t 

A whirling speed with respect to the- rotating coordinates thus corre- 

sponds to a whirling speed ccf = c% + cjq with respect to the fixed coor- 
dinates. A shaft critical vibration corresponds to = 0 in the 
rotating coordinate system or to oif = cd in the fixed coordinate system. 


Example of Rotor With Locked Hinges 


An example that involves a partial use of complex variables is given 
on page 255 of reference 2. The problem given there of a mass particle 
moving on the inner surface of a rotating spherical bowl is mathematically 
equivalent to the disturbed motion of a flywheel and shaft or of a rotor 
with locked hinges. The equations of motion obtained in real form in 
rotating coordinates 


*a " 2ay a - = 

y a + 2<m a - o£y a = 


B a x a K 
M “ M Xa 

► 

Bay a K 

M “ M y a 


( 1 ) 


are combined in the single equation 
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where 


z a = x a + iy a 


is the complex position vector in the rotating coordinate system. The 
complete solution, if small damping is assumed, is 


= cie- v l t+1 “^ t + c^Z*-**** (3) 

where 



The path of the motion is represented by rotations of a complex vector 
in a plane. 

The use of a complex variable has thus cut in half the number of 
equations to be handled and has yielded a solution from which the geo- 
metric path of the motion may easily be reconstructed. The advantage of 
the z-notation is not fully realized, however, unless it is used from 
the very beginning of the problem. The close similarity of this problem 
to the rotor-vibration problem makes it worthwhile to show the full appli- 
cation of the z-notation to the preceding example. The complex variable 
z a at any instant determines the position of the mass particle relative 
to the rotating coordinate system. If the position in a fixed coordinate 
system is denoted by zp. 


z f 


z a e iajt 


( 4 ) 


and z a can be treated as a generalized coordinate in the Lagrangian 
equations of motion. The kinetic- and potential-energy expressions can 
be immediately written as 


T = | Mz f z f 



+ icoz 



a " 
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V = \ Kz a z a 

A dissipation function for damping that depends upon motion relative to 
the rotating system can be written 


The equations of motion are now obtained by considering z a and z a as 
generalized coordinates in the Lagrangian equations. Substitution in the 
equation 


a /a t \ _ 8t_ + 8f_ + 8v_ _ o 
dt ^ a / d5 a al a Sz a 

thus yields the equation previously given 



The same method can be used to obtain the equations of motion in 
the fixed coordinate system. In this case, 


T = i MzpZf 


V = | KZpZf 


F = 2 B a ( if “ icoz f) ( z f + i“ z f) 
The equation of motion in terms of Zf becomes 


( 5 ) 


z'f + y ( £f " iu)Zf ) + | z f = 0 


( 6 ) 
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and the solution for small values of damping is 


z f 


C^e 




t+i^K/Mt. 


+ C2e 




,t-i)p0it 


( 7 ) 


This solution shows that the motion consists of two circular vibrations 

in opposite directions and, moreover, that for U3 > \j K/M the first term 
represents unstable motion > that is ^ the vibration has negative damping. 


This example illustrates a shaft critical speed for m = ^K/M and 

a self -excited instability for 03 > A discussion of the physical 

picture of this instability due to damping is given on page 293 of 
reference 1 . 


The effect of damping in a nonrotating part of the system can be 
included in the analysis merely by adding to the previous dissipation 
function the term 


2 BfMf 


The equation of motion then becomes 


MZf + B f z f + B a (zf - io3z f ) + KZf = 0 


( 8 ) 


The solution for small values of damping becomes 


z f = C ie 


Bf Eg/ 

2M 2M 


03 


^k/mj 


1+i 


i^K/M 


+ C2e 


2M 2M 


CD 


{k/m 




( 9 ) 


The motion is now unstable above the speed 


Bf\ 


03 = \|£|1 + -pi 

B a/ 


( 10 ) 


Hinged Rotor 

Inclusion of the effect of hinge motion in the plane of rotation 
increases the number of degrees of freedom and the number of equations 
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of motion. For example, three hinged blades and two directions of pylon 
deflection give five degrees of freedom to be considered. If special 
linear combinations of the hinge deflections 0k are used as generalized 
coordinates, no more than four degrees of freedom need be considered 
simultaneously. The use of complex variables reduces these four equa- 
tions to two equations . 

Appropriate variables in the rotating system for a three -blade 
rotor are 


e o 

= 7' 

(p 0 + 

0i + 02 ) 



i 

f 

2rti 


9 1 

II 

V p o + 

0]_ e 5 + 

0 2 e 


i 

( 

4rti 


0 2 


fo + 

0qe 5 + 

02 e 


4rti' 
5 ] 


Ml 

5 I 


(ID 


These variables and their complex conjugates satisfy. the relations 

§1 = -e 2 
0 2 = - 0 ! 


and also 


0q9q + ^2^2 ~ + 

= ^"( p 0 2 + Pi 2 + P 2 2 ) 

The variables 0k > by virtue of their meaning, are referred to a rotating 
coordinate system. The special linear combinations of the 0k denoted 
by 0^ are also referred to a rotating coordinate system. The appro- 
priate variables to represent the hinge deflections when fixed coordi- 
nates are used are defined by 


t k = e^ 1 "* 


and ^ is the complex conjugate of 


( 12 ) 
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Geometrically, 0^ or ^ is the complex vector representing the 
displacement due to hinge deflection of the center of mass of all the 
blades, just as z represents the position of the shaft due to pylon 
deflection. It will be shown later that, in the equations of motion, 

©1 is coupled with z and 0 q is an independent principal coordinate . 

Equations (ll) when solved for (3q, Pi, and p 2 become 


9 0 + 9 1 + 0 2 = bi Po 


0 o + 0ie" ia + 0 2 e ia = biPi 
0 q + 0j_e^- a + 02e“^ a = biP2 
Then, in a mode involving 0 q_, 


0 O = 0 

0 1 = | e 1 ^ 

0 2 = -Si 

P 0 = sin (Ogt 
Pi = sin (coat - a) 
P 2 = sin(aoat + a) 


(15) 


Equations (13) show that, in the 0i-mode, the blades are undergoing 
sinusoidal vibrations 120° out of phase with one another in a manner 
analogous to three-phase electrical currents. 
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General formulas for any number of blades are 


a _ bi 

a i — 
«J n 


n-1 


Pk e 


ijak 


k=0 


a = 


2* 

n 


9 n-j 


0 n ~ 9 0 


n-1 


k=0 


o n-1 

•A -^51 si 2 


n 


k=0 


( 14 ) 


Derivation of Equations of Motion 

The equations of motion and the characteristic equation of whirling 
speeds are herein derived for the general case of three or more equal 
blades on a pylon that may have different stiffness properties in dif- 
ferent directions of deflection. The effects of damping in the blade 
hinges and in the pylon are included. The equations are first formulated 
in a nonrotating reference system. The required modifications are then 
given for the case of isotropic support stiffness. The corresponding 
equations referred to the rotating coordinates are then obtained. 

Let the position of the center of mass of the kth blade be repre- 
sented by the complex quantity z^- in the plane of rotation. (See 
fig. 1.) Let the bending deflection of the pylon be represented by zf 
in a nonrotating coordinate system and let (3^ be the hinge deflection 
of the kth blade. Then 


z k = z f + (a + 
The complex velocity is 


Zk = Zf + jbi(3 k e i ^k ; + i^a + be'*'^ 


>i (ak+cut ) 


( 15 ) 


(16) 
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Because only small displacements are being considered, the exponential 
factors containing 0^ can be expanded and only the terms that lead to 
quadratic terms for the kinetic -energy expression need be considered. 

Some terms can be ignored either because they cancel after summation 
for al l the blades or because the corresponding derivative expressions in 
the Lagrangian equations vanish. The substitution 


e iPk 


= 1 + iPk - 



leads to the following expression for the kinetic energy of translational 
motion of the kth blade: 


Ik - 2 ^^k-^k 


(IT) 


where 


£ k t k = i f i t + 5 f bi(p k + + 

Zf(-bi)^0 k - icu0 k ^e"^ a ^ +aJ ^ + b^p k ^ - co^abp k ^ 

The kinetic energy of rotation about the center of mass of the blade is 


T r = | V 2 0 k 2 


(18) 


The effective mass of the pylon may be different in the Xf- and in the 
yf- directions . Allowance for this possibility is made by writing the 
kinetic energy of the pylon as 


T s = 



Vf 


mfifZf + Arrif 


Zf 2 + Zf 2 | 


( 19 ) 
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where 


m-f 


m X + Uly 
2 


Anif = 


m x — it^- 

2 ~ 


The total kinetic energy is the sum of the expressions for the separate 
kinetic energies . 


The pylon spring constant may differ in the Xf- and in the 
yf-directions and, consequently, the potential energy can be expressed 
as 


V = -iKfZfZ^ + AKf 


2 - 2 
Zf + Zf 

2 




( 20 ) 


The effect of damping will be expressed with the aid of a dissipa- 
tion function. If damping exists in the pylon, in the rotating shaft, 
and in the hinges, this function becomes 


F = | (hf Zf if + ABf Zf 2 ± Zf2 + B a z a z a + XI B p3 k 2 ' 
' d k=0 


( 21 ) 


The sum of the various energy expressions for all the blades, 
expressed in terms of the variables Zf and in the nonrotating 

coordinates, becomes 
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The Lagrangian equations of motion are 


d /dT \ 
dtydzf J 



8t i 8f ! 8v _ q 

dzf 8zf 8zf 


> 



(23) 


and similar expressions for the other variables. The equations of motion 
in fixed coordinates then become 


(n + nmb)z’f + Bzp + B a ^Zf - ia>z^ + Kz^ + Amz-p + ABzp + 
AKz f + nm^-L = 0 


nn^Zj + 2nm b 


1 + hlv i ■ 2ia * 


ntbb 


K, 


(3 


= 0 


&1 ~ °^l) + ^(&1 " ia *l) 


nnib 


2 a » 

^ v + 


y 1 + jl)(^ - 21a 4 - <^k) + ^s(4 - la <k) + 


K P 


b ik T ^ ?k 


= O' 


(24) 


where ^ refers to the ^-variables other than ^ and £ n _i* The com- 
plex conjugates of these equations are also obtained but give no additional 
information. Each complex equation is, of course, equivalent to two real 
equations. It is noticed that the first two equations contain only the 
variables zp , zp, and £1 and that the third equation represents n - 2 

equations, each containing one independent principal coordinate The 

physical meaning of this partial separation of variables is that a blade 
motion represented by £q involves a motion of the common center of mass 
of the blades and, thus, a coupling effect with the pylon. Blade motions 
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in which the common center of mass does not move are represented by 
£q, ••• 5 n * For ' t ^ Lree 'blades, the only such mode is the one corre- 
sponding to £q. In this mode, all the blades move in phase; the motion 
is always damped and does not lead to instability. 

The equations of motion of a one- or two-blade rotor are somewhat 
different from equations (24) . The difference is connected with the 
circumstance that a rotor of three or more equal blades has no preferred 
direction in its plane; whereas, a one- or two-blade rotor has different 
dynamic properties in directions along and normal to' the blades. Only a 
brief statement and the final equation for shaft critical speed will be 
given for the one- or two-blade rotor. 

The equations of motion involving Zf and can be written more 

compactly by use of the notation 

D = — D 2 = 

dt dt 2 


and the substitutions 


5 £ = Xf 

M “ 


Be 




2\ = A P 


Bo 


= h. 


bit 


rs 

D j 


- AXj. 
M 1 


Kc 


mbb2 ( i + 3 


= a 2 


nmb 


+ nm^ 


= M- 
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Then 


D 2 + AfD + Ag^D 
hD 2 Cl = 0 


ia)) + i\ zf + (tt 1)2 + ^ + lr) sf + 


2 r^ 


— — — D 2 z f + KD - icu ) 2 + Ap(D - ito) + 

r 2\ L 


V (25) 


ao^A^ + Ag Cl — 0 


or, briefly, 

A^(D)zf + AA]_^(D)zf + Aj_2(®)Ci = 0 
A pi (D)zf + Agg(D)Ci = 0 


(26) 


The Characteristic Equation 


The general form of solution of equations (26) is an elliptic 
whirling motion that can be represented by 


Zf = + C 2 a- l3ift 


io>ft 


i f = + U 2 e 


Ci - + 046^ 


(27) 


Special cases of this motion include whirling in a circle (Cg = C 4 = 0) 
and linear vibration (Ci = Cg, Cj = C 4 ) . Substitution of equations (27) 
in equation ( 26 ) gives 
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jAu(icof)Ci + M 11 (i(0 f )c 2 + + Ja 1X (-io>f) 

Mu (-iflif ) C]_ + A i2 C4] e _i ^ >t = 0 

|a 2 i (iciif ) c i + A 22 (io^ Cjj e la>ft + |a 2 i (-i%) C 2 + 

A 22 (“ 1 *%) C4] e“ i ^' t = 0 


C 2 +| 


(28) 


In order for equations (27) to be a solution of equations (26), equa- 
tions (28) must be satisfied for each value of t. The coefficient of 

each time factor e^ 0 ^ or e must therefore separately vanish. 

Because each bracketed expression represents a complex quantity that 
vanishes, its complex conjugate also must vanish. The condition for a 
solution can therefore be expressed by the vanishing of the first 
bracketed terms and the complex conjugates of the second bracketed 
terms. Hence, 


A llW c l-i M n( i( 8f) c 2 + Ai 2 (io>f)c 5 = 0 

A2i(io)f)Ci + A 22 (ioif) C3 = 0 

Mu Cl + An (i^) C 2 + A12 (icof) = 0 

A 2 1 ^icn^ j C 2 + A 22 = 0 


(29) 


where A-j^im^ is the complex conjugate of A-q (-iffif ) and is obtained 
from A^_]_ (i<^f) by changing iu> to -io> without changing io>f. The 
characteristic equation giving the rotational speeds is the determinant 
of the coefficients of Ci, C 2 , Cj, and C4 equated to zero. With the 
second and third columns interchanged for symmetry, the determinant becomes 


All 

A i2 

a 21 

a 2 2 

Mu 

0 

0 

0 


Mn 

0 

0 

0 

All 

Al2 

^21 

a 22 


( 30 ) 
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The expanded form of this determinant is 


where 


( A 11 A 22 “ A 12 A 2l) ( A 11 A 22 " ^12 A 2l) “ ^11 AA 11 A 22 A 22 “ 0 (51) 


A n = -% 2 + iai f A f + 1A a(“f - ®) + ^ 

An = “^f 2 + im f*f + ~ 

— Amf o AKf 

AA-n = M 1X = - -p- oif^ + icof AXf + 


A 12 A 21 = a 12 a 21 = ' T ^V.2\ ^ = *3^ 


2 1 + 


a 22 = “(^f - + iAp(tt>f - cd) + cd^Ai + A 2 

X pp = -(clip + 0>) 2 + iAp (oif + 03^ + 03 2 Ai + A 2 


The roots oif of this equation are the characteristic whirling speeds 
of the rotor. 

For the case of isotropic supports, 


Mu = 0 

and the equations of motion are satisfied by equations (27) with 
C2 = C4 = 0. 

The characteristic equation is then simply 


a 11 a 22 ~ A 12 A 21 " 0 


( 32 ) 
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In a rotating coordinate system, the complex coordinates are z a and 
where 

H - V 1 ”* 

El - 

Then 

Dzj = ^Dz a + iti)z 
D£i = (d9 x + imeje 10 ^ 

If the whirling speed in rotating coordinates is represented by <u a , 

z a = 

9i . 

The characteristic equation is then obtained by substituting + <d 
for : 

All(c% + <»)a 22 (c% + o>) - A 12 (“a + “) A 2l(<% + <*>) = 0 ( 33 ) 

The characteristic equation can thus be stated in terms of a whirling 
speed in either the fixed or the rotating coordinate system. 

METHOD OF APPLYING THEORY 


Application Neglecting Damping 

In plotting curves for use in applications of the theory, it is con- 
venient to consider one of the pylon bending frequencies = \j K X ^M X as 
a reference frequency and to refer all other frequencies as well as the 
rotational speed co to the reference frequency as unit. The number of 
independent parameters is thus reduced by 1. All quantities in equa- 
tions (31) to (33) are then expressed nondimens ionally. 
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The natural whirling speeds and the three types of vibration - ordi- 
nary, self excited, and shaft critical - can now be predicted from a study 
of the roots of equation ( 31 ) in which oif is considered a function of u> 

for fixed values of the other parameters . 

The case of no damping will be considered first. Because equa- 
tion ( 31 ) with damping terms omitted is of the fourth degree in oif 2 

and of only the second degree in cjd 2 , it may be solved conveniently by 
first choosing values of o>f and then solving the equation for u> 2 . 
Similar indirect methods can be used with equations (32) and (33)- 
Special methods to be used when damping is included will be discussed 
later . 

The meaning of equations ( 31 ) to (33) will be illustrated by exam- 
ples . The real part of o>f. will be plotted against o> for selected 

values of the parameters A;p, A2> A3, and s. The simplest case is 

that in which the mass of the blades is so small that any force on the 
pylon due to blade motions is negligible. The pylon motions are then 
independent of the blade motions. This case is obtained by putting 
A 3 = 0. The characteristic equation ( 31 ), (32), or ( 33 ) then factors 

into expressions yielding straight lines and hyperbolas. 

An example of a rotor with particular values of the parameters is 
plotted as long -dash lines in figure 2. The horizontal straight lines 
correspond to pylon bending and the slanting hyperbolas correspond to 
hinge deflection. Each curve represents the trend of one of the real 
roots o>f. As A 3 increases slightly from zero, the greatest changes 
in the curves occur in the vicinity of the intersections of the straight 
lines with the hyperbolas. Here each branch breaks away from the inter- 
section and rejoins the other branch. At a gap, such as C in figure 2, 
the number of real roots of the frequency equation is reduced by 2. The 
missing roots are complex conjugate numbers, and one of them must have a 
negative imaginary part, which implies a self -excited vibration. 

Consider the interpretation of figure 2 as m is gradually increased 
from zero . At zero rotational speed, the values of cpp are the natural 
frequencies that could be excited as ordinary vibration by applied vibra- 
ting force. Positive and negative . values occur in pairs of equal magni- 
tude and correspond to linear vibration modes represented in complex nota- 
tion as 

Zf = + e-^) 

As cd increases from zero, the positive and negative values of o>f 
no longer are equal in magnitude. The normal modes are therefore whirling 
motions with angular velocities equal to the plotted values of o>f. 
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The shaft critical speed is the rotational speed at which o)f = a) 
and hence is given by the point A where a 45° line through the origin 
intersects the a>f-curve. This speed corresponds to the peak for vibra- 
tions excited by imbalance in the rotating system. As a) increases 
above the shaft critical speed, the modes of whirling are stable until, 
for the case of no damping, the value of <Df becomes complex at the 

value of a) at which a vertical line is tangent to the plotted curve . 

This point B is the beginning of the self -excited range. At the point D, 
the motion again becomes stable. The real part of o>f is plotted ir the 
region C as a short-dash line. The complex roots in the region C have 
been calculated and plotted in figure 3* 

The point E, at which oif = 0, is of some interest. At this speed, 
a vibration of the blades could be excited by a steady force (u>f = 0, 
cua = -cd), such as the force of gravity if the plane of the rotor is not 
horizontal. 

Because the most important information to be obtained from the fre- 
quency equation is the critical value of 0) for the shaft critical and 
self -excited vibrations, a set of charts that gives this information for 
a large variety of values of the physical parameters has been prepared. 
These charts are given in figures 4 to 6, which correspond to values of 
stiffness ratio K y /K x = s of 1, and 0, respectively. The use of the 

charts is illustrated by a numerical example. Suppose the values of the 
parameters for a certain rotor are Ai = 0.07, A2 = 0.22., A3 = 0.1, 

s = 1, and aij- = 155 cycles per minute. A straight line, such as AB in 
figure 4, is first drawn to represent the function cd2a^ + A 2- This line 
intersects contours A3 =0.1 at 0)2=0.77 for the shaft critical point 

and a£ = 1.6 and 4.85 for the beginning and for the end of the self- 
excited range, respectively. With a reference frequency of 155 cycles 
per minute, these values correspond to actual rotational speeds of 136, 
196, and 342 rpm. 

All possible values of A]_, and A3 are thus covered by suit- 

ably changing the straight line AB. The general effect of the stiffness 
ratio s is not large; any case can therefore be estimated with a fair 
degree of accuracy by use of figures 4 to 6. 


Possibility of Avoiding Occurrence of Vibration 

Figures 4 to 6 can also be used for the inverse problem of finding 
the values of the parameters that are required to obtain given values of 
critical rotational speed. These figures show that to eliminate entirely 
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the self -excited instability requires that Aq be equal to or greater 
than 1 . The shaft critical instability can be entirely eliminated only 
with a value of Aq in a small range near 4 and with s = <*> or s = 0 . 
These values of Aq differ radically from present designs in which a 
typical value is 0.07* 

The satisfactory requirement of keeping the instabilities outside 
the operating range of rotational speed i s found by first picking a 
reasonable value of the pylon frequency ^K x /Mx to fix the scale unit 
for to and by then observing the combinations of Aq and Ap that 
can be used to avoid the critical Aj -contours . 


Effect of Damping 

The effect of damping has been included in equation ( 31 ) through 
the parameters Af, AAf , A a , and Ap. A method of computation similar 

to that used in flutter theory appears preferable to attempting to solve 

the equation directly for Of. The beginning and the end of an unstable 

range can be found by the following method: At a limit point between a 

stable and an unstable speed range, the value of cop is real. Equa- 
tion (31) is first separated into real and imaginary parts with Of con- 
sidered real. Each part is considered a functional relation between cof 
and 00 and is plotted for a given set of values of the parameters . The 
intersections of the real and the imaginary equations give the rotor 
speeds and frequencies corresponding to the beginning and the end of the 
unstable ranges. In the computations, it is preferable to choose values 
of o>f and to solve the equations for the corresponding values of to. 

The explicit form for computation in the simplest case of isotropic 
supports, with damping in the pylon and in the hinges but not in the 
rotating shaft ^Aa = 0^, is obtained from equation (32) rearranged as 
follows : 

For the real equation 


o)2 - 2 BrO) + Cr = 0 


( 54 ) 


CJD 


Br ± ^Br 2 - Cr 
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where 


Br = 


(&f 


1 - Ai 


1 + 


AfXp 


+ a) 


and 


C 0>f 2 [ 1 + Ag Ayaf 2 + Af Ap) 

1 " A 1 ( a * 2 _^2 + h . 

^ M 


For the imaginary equation 


cd 2 - 2Bjo> + Cj = 0 


(35) 


where 


and 


B I 


1 

1 - Ai 


Rf 


^ — f-03f 2 + 

2AfCDf \ 



C I = “ 


1 


1 - Ax 

A f 



2 + £c\ 

M/ 


+ (-(Up 2 



The real and imaginary equations for the most general case of equa- 
tion (31) can he written, respectively, as follows : 

A a 2 (l - Ai) 2 0) 6 + |Ri(l - A1) 2 + A^R^O) 4 + (RpRj - Iplj + A a 2 R 2 - R 5 )^ + 

R1R2 - Ipl 2 ~ R4 + Aj 2 ^® = 0 (36) 


Il(l - Ai) ; 


+ V*3 


03" 


+ . (RjIj + R3I1 + A a 2 I 2 - I^oj 2 + 


R^l2 + ^2*^1 “ I4. ~ 0 


(37) 
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where 


Rl ’ ifH^ 2 + Sjf 1 * 2 + 1|) " *m + ^X 7 * ^ + *») 


= &L .2 + 5lU>.5 + O + ^U 2 + &W»* & + 


M 


Mx, 


wv 


M 


M 


-CCf 11 




Pf A X 


M 


R2 = + A2^ * CDf^Xp^ 

I2 = 2 ^-oif + A2^GDf Aj3 

R^ = -2 ^1 - A]_^ ^ + A2^ - + Ap^ 

1^ = -2 ^1 - <x >£ Ap + 4o^»Ap 


R^ = 2oifTAj 


\ = 2o>f 4 A5 


j^Oif 2 + — ^ (-^f 2 + ^ 2 ) ” 9f ' 2 (v + \l) ^0 

^-aif 2 + — ^a^Ap + o>f (^f + A a ^ (-o>f 2 + A2^ 
-(! " Al) ^-oof 2 + - AaAp 


R5 = 2aUf H A^| 

I5 = 2tDf^Ajjj^l - A^aif^Af + A a ) - 2Aa0^j 


Examples of calculated cases with damping are shown in figures 7 
to 9. The presence of small amounts of damping in both the pylon and 
the hinge degrees of freedom does not greatly change the predictions 
that would be made from the equations with no damping. The plot of the 
real equation is practically the same as the plot obtained when damping 
is neglected. The intersections of the curves of the imaginary and the 
real equations with any reasonable value of Af/Ap are near the points 

that would be considered the limits of the unstable range if damping 
were neglected. Increasing the amount of damping decreases the gap 
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between the limits of stability until the unstable range is finally 
eliminated. An approximate solution for the amount of damping required 
to eliminate the self -excited instability is obtained by requiring that 
the damping be at least large enough to make the curve of the real equa- 
tion pass through the point where o>f = 1 and cd is the value given by 

the equation 


1 = 0 ) 


- \Ja£ Ai 


+ A2 


The values required in the case of s = » have been computed and plotted 
in figure 10. The elimination of self -excited vibration by damping thus 
looks promising and merits further study with reference to specific 
application . 


LIMITATIONS AND FURTHER DEVELOPMENTS OF THE THEORY 


Polar Symmetry 

An important idea in the rotor vibration theory is the concept of 
polar symmetry. This concept implies the absence of a preferred direc- 
tion in the plane of the rotor. A rotor of three or more equal blades 
has polar symmetry. A rotor of two blades or one with unequal centering 
springs does not have polar symmetry. A pylon for which K x = Ky, 

Bx = By, and m* = my has polar symmetry. The possibility of solving 

the rotor vibration problem in terms of exponential or trigonometric 
functions depends upon the existence of polar symmetry in the rotating 
parts or in the nonrotating parts or in both. The general case of no 
polar symmetry would lead to Mathieu functions or something similar. 


Two Blades 

A brief comparison between the two-blade and the general case is 
presented herein. Polar symmetry of the pylon is assumed. The shaft 
critical speed is obtained by substituting (% = 0 in the character- 
istic equation as expressed in a rotating coordinate system. For one 
or two blades, the equation obtained is 
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The first bracketed factor gives the beginning of a self -excited range 
and the second factor gives the end of the range. 

Equation (38) can be compared with the following equation for the 
shaft critical speed of a rotor with three or more equal blades and 
polar symmetry: 



A useful chart based on equation (39) Is given in figure 11; some experi- 
mental results of tests of a simple model are given in figure 12. These 
tests demonstrate the essential difference between the two-blade and the 
general case. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., August 24, 1956. 
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Figure 2.- The effect of coupling between pylon and hinge motions. 



Figure 3- - The complex frequency in the unstable range for = 0.07, 
A2 = 0.22, Aj = 0.1, and s = 1. 
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Figure- 7 .- Plot of real and imaginary equations for a typical case, s 
Ai = 0.07; A 2 = 0.22; and Aj = 0.198. 



Figure 8.- Plot of real and imaginary equations for case of s = 
A-l = 0.07, A 2 = 0.22, and Aj = 0.198. 
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CHAPTER II 

THEORY OF MECHANICAL OSCILLATIONS OF ROTORS 
WITH TWO HINGED BLADES 
By Arnold M. Feingold 


SUMMARY 


The mechanical stability of a rotor having two vertically hinged 
blades mounted upon symmetrical supports, that is, of equal stiffness 
and mass in all horizontal directions, is investigated arid reported 
herein. The frequency equation is derived and shows the existence, in 
general, of two ranges of rotational speeds at which instability occurs. 
The lower region of instability is bounded by two shaft critical speeds. 
At rotor speeds within this region, self -excited divergence of the rotor 
takes place analogous to the instability exhibited by a rotating shaft 
that is elliptical in cross section. Within the second instability 
range, the rotor system undergoes self -excited oscillations. Charts are 
presented giving the boundary points of both instability regions for a 
large variety of values of the physical parameters. The effect of 
damping is also included in the analysis. 


INTRODUCTION 


In reference 1, Coleman gives an analytical study of the mechani- 
cal stability of a rotor having three or more vertically hinged blades, 
mounted on flexible supports. It was shown that, in addition to the 
usual shaft critical speeds, self-excited vibrations occurred over a 
range of rotational speeds. Experiments with rotary -wing aircraft have 
confirmed the soundness of the analysis. 

The present paper is an investigation of the stability of the two- 
blade rotor mounted on symmetrical supports. As will be shown later, the 
results differ from those for a three -blade rotor. The reason for the 
different behavior lies in the inherent asymmetry of a rotor with only two 
blades. Motion of the center of mass of the blades of a two-blade rotor 
with respect to the rotor hub, due to small hinge deflections of the blades 
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must be normal to the line of the blades. This restraint, which does not 
appear in a rotor of three or more blades, results in the rotor system 
having different dynamic properties along and normal to the line of the 
blades. Therefore, with supports that have equal stiffness and mass in 
all directions attached to a two-blade rotor, two principal vibration 
axes of the rotor hub can still be distinguished. No preferred vibration 
axes can be distinguished for a three-blade rotor mounted on symmetrical 
supports . This distinction shows up physically in the shape of the vibra- 
tion modes. Whereas a three-blade rotor whirls in a circle, a two-blade 
rotor whirls in an ellipse, of which the principal axes are along and 
normal to the line of the rotor blades. 

A two-blade rotor can be expected to show, in addition to some fea- 
tures of a three-blade rotor, some of the characteristics of a rotating 
shaft that is elliptical in cross section. Such a shaft, mounted on 
symmetrical bearings, is known to have two critical speeds, which corre- 
spond to each of the two principal stiffnesses. (See, for example, ref. 2.) 
For all rotational speeds between the critical speeds, the shaft is unstable 
and diverges. It will be shown that an exactly similar phenomenon exists 
for a two-blade rotor. The existence of this region of instability for a 
two-blade rotor is predicted in reference 1, in which the formula for the 
shaft critical speeds bounding this instability range is given. In addi- 
tion to this region of instability, a second range of instability analogous 
to that exhibited by a three-blade rotor is also present. 

Only the case of symmetrical supports is analyzed in the present 
report. In the case of asymmetric supports, the equations of motion are 
linear differential equations that are difficult to solve because the coef- 
ficients vary periodically with the time (Mathieu type). Similar equations 
are obtained in the problem of a rotating elliptical shaft mounted on asym- 
metric bearings. (See ref. 2.) 


SYMBOLS 


a radial position of vertical hinge 

b distance from vertical hinge to center of mass of blade 

B damping force per unit velocity of rotor-hub displacement 

Bg damping force per unit angular velocity of blade displacement 

about hinge 

D time -derivative operator, d/dt 

F dissipation function 
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I moment of inertia of blade about hinge, m^b^ 

K spring constant of rotor-hub displacement 

Kp spring constant of blade self-centering spring 

m effective mass of pylon 

effective mass' of rotor blade 

M total effective mass of blades and pylon, m + 2ntb 

r radius of gyration of blade about its center of mass 

s arbitrary parameter 

t time 

T kinetic energy 

T 1^ T 2 kinetic energies of rotor blades 

T s kinetic energy of rotor hub 

V potential energy 

x,y displacements of rotor hub in rotating coordinate system 

X a ,Y a rotating coordinate axes 

Xf,Yf fixed coordinate axes 

x 0 ,y 0 values of x and y when t = 0 

x ±>yi displacement of first rotor blade in fixed coordinate system 
x 2 ,y 2 displacement of second rotor blade in fixed coordinate system 
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*(0]_ - P 2 ) 

2 



value of 9i when t = 0 


A = 


B 

Moip 



Al = 


hi + 


Kn 


*2 = 


Icd^ 


a 3 = 


“b 


M 1 + i- 

V h 2 y 


0) 


angular velocity of rotor (the dimensionless ratio ou/o^ is 
called oo in applications) 




03f 


003- 


natural frequency of rotor system observed in rotating coordi- 
nate system (used in nondimensional form in applications) 

natural frequency of rotor system in fixed coordinate system 
(nondimensional in applications) 


reference frequency, \f k/m 


MATHEMATICAL ANALYSIS 


Four degrees of freedom of the system are considered - horizontal 
deflection of the rotor hub in the x- and y-directions, and hinge 
deflections p^ and p 2 "the blades in the horizontal plane of the 

rotor hub. The rotor is assumed to rotate at a constant velocity a>. 
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Deflection of the rotor hub may be due either to the bending of a 
flexible pylon or to a rocking of the rotor craft upon its landing gear. 
Ground-resonance vibrations usually involve landing -gear flexibility. 

The mathematical treatment is the same in both cases, but the values of 
several of the physical parameters will depend upon which mode is being 
investigated. Throughout this paper, the terms "rotor supports" and 
"pylon" will be used interchangeably to denote the nonrotating structure 
coupled with the rotor blades . 

The mathematical treatment herein differs from that in reference 1, 
in which are used the complex notation and the notion of "whirling speeds," 
that is, directional frequencies resulting from the use of complex numbers. 
Although the method of reference 1 is valuable for systems, such as the 
three-blade rotor on symmetrical supports, which have circular modes of 
vibration, it offers little advantage for the present problem, in which 
the rotor performs elliptical motion. Rectangular coordinates accordingly 
are used in the present report and frequencies are used instead of whirling 
speeds. In comparing the results of the present report with those of ref- 
erence 1, care should be taken to distinguish between frequencies and 
whirling speeds. Whirling speeds have directional significance; whereas 
frequencies are essentially positive quantities and do not give any imme- 
diate information concerning the direction of whirl of the vibration. 

The equations of motion are set up in a coordinate system rotating 
at the velocity co. Let the deflection of the rotor hub be represented 
by x and y in rotating coordinates. (See fig. 1 in which the inter- 
section of the coordinate axes represents the undisturbed position of the 
rotor hub.) The disturbed positions of the two blades in fixed coordinates 
are 


X 1 = (x + a + b cos P]_)cos cut - (y + b sin Pq)sin cot 


yq = (y + b sin p^)cos cot + (x + a + b cos p 1 )sin cot 


and 


X 2 - (x - a - b cos P 2 )cos cot - (y - b sin p 2 )sin cot 
y 2 = (y - * s i n p 2 ) cos cot + (x - a - b cos p 2 )sin cot 
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The kinetic energies of the two rotor blades are 

T 1 = ^b^l 2 + ^l 2 + r2 (“ + Pi) 2 ] 

and 

t 2 = ^b^ 2 + ^2 2 + + P 2 ) 2 ] 

The kinetic energy of the pylon is 

T s = |mji 2 + y 2 + o£(x 2 + y 2 ) - 2co(xy - xy] 

Because only small displacements from the equilibrium position are 
considered, the trigonometric expressions containing (3^ and 3 2 may be 

expanded as power series and only the terms that lead to quadratic terms 
in the energy expressions need be retained. Thus 


and 


sin 3 X = 3 -l 

3 1 2 

cos 3q = 1 - — 


sin 3 2 = 3 2 

3 2 2 

cos 3 2 = 1 - 

New variables introduced to replace and (3^ are 


00 = |(Pl + P 2 ) 

01 = I(P1 - p 2 ) 


where 0-|_ represents the shift, due to hinge motion, of the center of mass 
of the two blades with respect to the rotor hub. The introduction of 0q 
and 0-|_ results in a partial decoupling of the equations of motion. 
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The total kinetic energy of the system is 


T = T]_ + T 2 + T s 


Only the quadratic terms will he retained in the kinetic -energy expression, 
because the terms of lower degree vanish in the Lagrange equations of 
motion. Then 


T = |m|x 2 + y 2 - + co 2 (x 2 + y 2 ) - 2cu(xy - xyTj + 2y0;L 


2m 2 y0^ + 


2osc0^ 


2 ux 03_ + 





+ 0 - 


I 


The potential energy of the system is 

V =. 1k(x 2 + y 2 ) + + p 2 2 ) 

■ §K(x 2 + y 2 ) + i(e 0 2 + 01 2 ) 


Two types of damping of the rotor system are assumed to exist: 

(l) damping in the rotor supports, which is proportional to velocity 
displacements of the rotor hub in a fixed coordinate system, and (2) 
damping in the blade hinges. The dissipation function F then becomes 



, 2/^2 


or(x' 


+ 


y 2 ) 


2u)(xy 


xy 


3 + 


’0 



where B is the damping force per unit velocity of rotor-hub displacement 
and Bp is the damping force per unit velocity of a rotor-blade displace- 
ment about the blade hinge . 
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The Lagrange equation of motion for the variable x is 


_d/dT\ _8 T, 8 v. 8F = n 
dt \8x / 8x dx dx 

The use of similar expressions for y, 0 O , and 0-|_ lead to the following 
equations of motion: 


(d 2 - m 2 + — D + — ^ 

V m m; 


x TT™ D^l " D + rru>)y = 0 

M x 


( 2a) D + |^)y = 


(1) 


2m Dx + 


[i + ^rP + Jg- D + J5L + ^ 


b 


0q + (d 2 - m 2 )y = 0 (2) 


( 2cu D + |“) x + ^( d2 " uj2 ) 0 l + ^D 2 - m 2 + | D + J^y = 0 (3) 


'(* • 


|2 + D + a? | 

mbb 2 n^b 2 


0 O = ° 


(4) 


where 


D = — 
dt 


1)2 = JL 

dt 2 


Equation (4) can be solved independently of the others because it 
is an equation in only one variable 0 q . Equation (4), which also was 
obtained in the study of the three 4b lade rotor (ref. 1 ), represents blade 
motion with the blades moving in phase, uncoupled with pylon motion. 
Motion in this mode is damped and does not lead to instability. 
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Assuming solutions of the form 


x = x Q e 




y = y Q e 


la ^ 


9 1 = 9 1 e 




( 5 ) 


and substituting these solutions in equations (l) to (3) gives the charac- 
teristic or frequency equation 


-tc^ 2 -tu 2 + i Aj^+ 1 4 a 


3 °*% 


2oiJ^-iAco 


2o ix% -a) a 2 +iApa^+A2+A]_0D 2 -av^-cu 2 


2axD a -iAa3 -u^ 2 -a^+iAa^+l 


= 0 


( 6 ) 


where the nondimensional parameters 


Ai = 


b 1 + 


Ap = 


Kr 


Ia^ 


a 3 = 


“b 


Mil + — 

\ t^2 , 
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A = -S- 

M03j. 


A p -A. 

P IOJp 


have been introduced, and the rotational velocity <n and t he frequency 
have also been made nondimens ional by using = ^75 as reference 
frequency . 

The frequency generally is a complex number, of which the real 

part is the frequency of the vibration and the imaginary part determines 
the rate of damping of the vibration. If the imaginary part of is 

negative, the vibration increases in amplitude with time and the rotor is 
unstable. 


DISCUSSION OF FREQUENCY EQUATION 


Case of Zero Damping 

If the damping parameters A and A^ are neglected, the frequency 
equation (6) may be expanded to 

2A^(a£ + t% 2 - lj + po £<% 2 - + 

c% 2 - l) 2_ l l-2A 5 ^a> 2 + <% 2 + l) + c% 2 - A 2 - Aqtopj = 0 


(7) 


where is the natural frequency of the rotor system in a coordinate 

system rotating with the rotor. (Although equation (7) is a cubic 
equation in both o? and , rectangular hyperbolas of the form 
0^2 = o? + s, where s is an arbitrary parameter, intersect equation ( 7 ) 
at only two values of . For purposes of computation, therefore, equa- 
tion ( 7 ) can be reduced to a quadratic equation in o£ by replacing 
with o£ + s . ) 
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The solutions for zero damping (equation (5)) represent motion of 
the pylon in an ellipse expressed relative to the rotating coordinate axes. 
In fixed coordinates , the pylon would move in an ellipse precessing at the 
velocity oo. This motion can be resolved into simultaneous circular motion 
at the two frequencies | cjo + cnj and |o> - ofc|, in which the vertical 

lines indicate that the quantity inside is to be considered positive. If 
the pylon is subjected to a harmonic force in the fixed coordinate system 
of frequency aif, resonance will occur at each of the frequencies 


o>f = |o) ± o^l 

The frequency o>f will be referred to as the natural frequency of the 
rotor system in fixed coordinates. 

The graph of the frequency equation (j) for a typical set of values 
of the parameters is given in figure 2 in rotating coordinates and in 
figure 3 in fixed coordinates. For zero coupling between the blades and 
rotor hub - that is, when Aj equals zero - equation (7) factors into 

straight lines and a hyperbola, which are shown as long-dash lines in 
figures 2 and 3. The straight lines represent hub motion and the hyper- 
bola represents blade motion. A small increase in results in a 

breaking away of the curves at their intersections to form two self- 
excited regions. It is interesting to compare figure 3 with figure 4, 
which is the graph of the natural frequencies of a three -blade rotor 
having the same values of Aj_, Aq> and Ay 

The shaft critical speeds, or natural frequencies that would be in 
resonance with an unbalance in the rotor system, are found by putting 
= 0 in equation (7)* Figure 2 shows two such speeds, at points A 

and B (shown also in fig. 3)> that bound a region in which is a 

pure imaginary number. If cc^ is a complex root of the characteristic 
equation, the complex conjugate of will also be a root and one of 

the two roots will have a negative imaginary part implying instability. 
The rotor system will thus be unstable for all rotational speeds between 
the two shaft critical speeds. Because is a pure imaginary number 

in this region, the frequency of the resultant self -excited vibration is 
zero in a rotating coordinate system - similar to the shaft critical 
speeds - and will appear as a self -excited divergence of the rotor. 

The equation of the shaft critical speeds is 

|(l " “ 2 ) ( A 2 + A l 0 ^) - 2^] (1 - <£) = 0 


( 8 ) 
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The first factor gives the lower shaft critical speed. The second factor, 
which depends on only the reference frequency, marks the end of the range 
of instability and is the second shaft critical speed. Formula (8) and 
an experimental verification of it are given in reference 1. A convenient 
graph of equation (8) is given in figure 5. It will be noticed that it is 
impossible to remove the two shaft critical speeds or the instability 
region between them by any possible change of the parameters Aq, Ag , or 

Ay, that is, without the introduction of damping, self -excited vibrations 
will always occur below the rotational speed o^. 

Instability also occurs in a range of rotational speeds greater than 
(Op. This range is shown in figures 2 and 5 as the region bounded by the 

points C and D and is similar in origin to the self -excited region 
exhibited by the three-blade rotor. In this region, the roots of the fre- 
quency equation are complex and self -excited vibrations will take place. 
Unlike the three -blade rotor, however, the rotor hub will be seen from a 
stationary position to be simultaneously executing self -excited vibrations 
at two different frequencies. Physically, of course, the rotor is moving 
in an ellipse at the frequency while precessing at the velocity cu. 


A chart showing the lower and upper limits of this instability 
region for a wide choice of values for the parameters Aq, Aq, and Aj 

is given in figure 6. The chart is used by drawing a straight line that 
represents the function (l - k-A^a? plotted against Aqo£ + Aq • The 
intersections of this straight line with the proper Aj-curves give the 

desired values of <x>. The short-dash line on the chart illustrates the 
method for the parameters of figure 2. 


The position of the instability region is very sensitive to the 
values of Aj. (See fig. 7-) As Aj increases, the region of insta- 
bility occurs at greater rotational speeds and moves to infinity for 
A 3 = jkl - Aq). For values of A^ greater than i(l - Aq) - that is, 


when the total effective mass of the rotor blades is greater than the 
effective mass of the rotor supports - the self -excited region does not 
appear . 


At certain rotational speeds, cpf = 0. At such speeds resonance may 

be excited by a steady force, constant in direction, acting on the pylon 
or blades - for example, gravity acting on a tilted rotor. The two-blade 
rotor has two such speeds, shown as points E and F in figures 2 and 3 . 

The mathematical condition for such points is that = 0 or = a£ 

in equation ( 7 ). The equation giving the rotational speeds at which this 
condition may occur is 
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Ag j {ko£ - lj - a? Ja^ ^ 



= 0 


(9) 


Equation (9) is plotted in figure 8, which is used similarly to figure 6. 
If A 2 = 0, then point E occurs at cu = 0. If A-, ^ i(l - At ) , as 

5 " V 1 

pointed out in reference 3, then point F occurs at oo = <*>. 


Effect of Damping 

The effect of damping will he determined in the same manner as for 
the three -blade rotor in reference 1. When the damping parameters A and 
Ap are retained^ equation (6) in expanded form can be separated into 

powers of having real and imaginary coefficients . ■ The terns of equa- 

tion (6) with real coefficients are 


2A^(o? + cD a 2 - l) + - 

(a) 2 + C% 2 - l) 2 J ^2A 5 (cd 2 + (% 2 + 

A 2 (-OXj 2 + Ag + AyO) 2 ^ (u£ - U ^ 2 ^ - 


l) + °h. 2 - Ag - Aja^j + 

2 AAp<% 2 (o) 2 - o^ 2 + l) (10) 


The terms with 


imaginary coefficients are 


i7 W 


2 ^ [P - + X ) ( A l^ 2 + *2 ~ % 2 ) 


V 

a 3 (3^ - 2o) 2 o^ 2 - o^jj + + cu^ + 1 - 

2 o) 2 o ^ 2 - 2 o ^ 2 - 2 o ? - A 2 ^ 2 - 


+ 


( 11 ) 
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At a boundary between stability and instability, is real. 

Such points are the intersections of the equations formed by setting 
expressions (10) and (11) separately equal to zero and plotting them on 
the same coordinate axes. Figure 9 shows a calculated case of damping. 
The imaginary equation is plotted for several values of the ratio of the 

damping parameters A^Ap, with A 2 assumed to be negligible. It is 
seen that, for large values of A^Ap, the boundaries of the higher range 

of instability are not far different from the boundaries found by neg-. 
lecting damping. Small values of A^Ap - that is, when most of the 

damping is concentrated in the blade hinges - lead to a beginning of the 
instability at lower rotational speeds. 

For small amounts of damping, the plot of the real equation is prac- 
tically the same as when damping is neglected. By introducing sufficient 
damping, however, the higher instability region may be eliminated. (See 
fig. 10.) The two shaft critical speeds and the instability region 
between them can also be removed by putting enough damping into the rotor 
supports, although a large amount of damping is required. 


Brief Description of Vibration Modes 

If damping is neglected, the shape of the free vibration modes can 
be found from the equations of motion (equations (l) to (4)) and the 
form of the solution (equation (5)). The rotor hub generally moves in 
an elliptical path in rotating coordinates although, at certain speeds, 
the motion may become circular or linear. At zero rotational speed, two 
of the three modes involve hub motion normal to the line of the blades , 
with concomitant blade motion. In the third mode, the blades do not 
move about their hinges and the rotor hub moves in a straight line parallel 
to the line of the blades at a frequency equal to . 

At the first shaft critical speed, the rotor hub diverges in a direc- 
tion normal to the line of the blades; whereas, at the second shaft critical 
speed, the hub diverges parallel to the blades . 

The forced responses of the system to a vibrator attached to the 
pylon can also easily be determined and show that those responses lying 
closest to the lines cop 2 = 1 are the strongest. When the coupling param- 
eter Aj is zero, no response occurs along the lines uif = | 2cjd 1 l|. This 

last conclusion is, of course, necessary if the theory is to give the cor- 
rect results for the degenerate case of massless rotor blades . 
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CONCLUSIONS 


The mechanical stability of a rotor having two vertically hinged 
blades mounted upon symmetrical supports has been investigated and 
reported herein. This investigation indicated that the main features of 
such a rotor system may be summarized as follows : 

1 . The vibration modes are generally elliptical, as opposed to cir- 
cular for the three-blade rotor. The ellipse precesses at a speed 00 
as observed from a fixed position; the result is six reso nan t or natural 
frequencies in a fixed coordinate system for a given rotor speed as 
against four natural frequencies for the three-blade rotor. 

2 . The asymmetry of the two -blade -rotor system gives rise to a range 
of rotor speeds in which self -excited divergence of the rotor occurs . 

This instability region is bounded by two shaft critical speeds. A 
three-blade rotor, in contrast, has only one shaft critical speed with 
no associated instability region. 

3. The two-blade rotor has a second region of rotational speeds at 
which self -excited vibrations occur. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., July 27 , 1956. 
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Figure 2.- Natural frequencies of a two-blade rotor in rotating 
coordinates for case of A]_ = 0.05, Ao = 0.20, and Ax = 0. 
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Figure 3 .- Natural frequencies of a two-blade rotor in fixed coordinates 
for case of A x = 0 . 05 , A% ~ 0 . 20 , and * ~ 
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Natural frequencies of a three -blade rotor in fixed coordinates 
for case of A x = 0.05, Ag = 0.20, and A, = 0.10. 
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Figure 6.- Stability chart for second instability region 
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CHAPTER III 

THEORY OF GROUND VIBRATIONS OF A TWO-BLADE HELICOPTER 
ROTOR ON ANISOTROPIC FLEXIBLE SUPPORTS 
(REVISED) 

By Robert P. Coleman and Arnold M. Fei'ngold 
SUMMARY 


An extension of previous work on the theory of self -excited mechanical 
osci llat ions of hinged rotor blades has been made. Previously published 
papers cover the cases of three or more rotor blad.es on elastic supports 
(such as landing gear) having either equal or unequal support stiffness 
in different directions and the case of one- or two-blade rotors on sup- 
ports having equal stiffness in all horizontal directions. The missing 
case of one or two blades on unequal supports has been treated. 

The mathematical treatment of this case is considerably more compli- 
cated than the other cases because of the occurrence of differential equa- 
tions with periodic coefficients. The characteristic frequencies are 
obtained from an inf inite -order determinant. Recurrence relations and 
convergence factors are used in finding the roots of the infinite 
determinant . 

The results show the existence of ranges of rotational speed at 
instability occurs (changed somewhat in position and extent) similar 
to those possessed by the two-blade rotor on equal supports. In addition, 
the existence of an infinite number of instability ranges which occurred 
at low rotor speeds and which did not occur in the cases previously treated 
is shown. 

Simplifications occur in the analysis for the special cases of infi- 
nite and zero stiffness in one of the axes. The case of infinite stiffness 
in one axis is also of special interest because it is mathematically equiv- 
alent to a counterrotating rotor system. A design chart for finding the 
position of the principal self -excited-instability range for the case of 
infinite support stiffness in one direction is included for the convenience 
of designers. It is expected that designers will be able to obtain suf- 
ficiently accurate information by considering only the cases of infinite 
and zero support stiffness along one direction together with the cases 
treated previously. 
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INTRODU CT ION 


It is known that rotating -wing aircraft may experience violent vibra- 
tions while the rotor is turning and. the aircraft is on the ground. It 
has been found that these vibrations .can be explained without considering 
aerodynamic effects and that they are due to mechanical coupling between 
horizontal hub displacements and blade oscillations in the plane of rota- 
tion. A theoretical analysis of this vibration problem is given in ref- 
erences 1 and 2. Reference 1 deals with rotors having three or more equal 
blades on general supports and reference 2 deals with two-blade rotors on 
supports having the same stiffness in all directions. 

Although in actual two-blade rotary-wing aircraft, the stiffness 
of the supports along the longitudinal direction is certainly different 
from the lateral stiffness, the equality of the stiffnesses was assumed 
in reference 2 because it permitted the mathematical simplification of 
dealing with differential equations having constant coefficients and it 
was believed that a theory employing such an assumption would be suffi- 
cient to indicate the nature of the most violent types of ground 
instability. 

The present paper gives a theoretical investigation of the general 
case of a two-blade rotor mounted upon supports of unequal stiffness 
along the two stationary principal axes . It thus generalizes the problem 
of reference 2, and rounds out the studies of ground resonance begun in 
reference 1. As was shown in reference 2, a two-blade rotor possesses 
cLifferent dynamic properties along and normal to the line of the blades. 
Equations of motion with constant coefficients for the problem treated 
in reference 2 could be obtained by using a coordinate system rotating 
with the rotor. This procedure succeeded because the supports were 
assumed isotropic (equal stiffness in all directions) . When the supports 
are anisotropic, however, it is impossible to avoid the appearance of 
periodic coefficients in the equations of motion. 

The present method of solving the differential equations of .motion 
follows closely the process employed in reference 5 for a vibration prob- 
lem in two degrees of freedom. The form of solution is expressed by an 
exponential factor times a complex Fourier series. Substitution of the 
formal solution into the equations of motion yields an infinite set of 
algebraic equations and an infinite -order determinant for the determina- 
tion of the Fourier coefficients and the characteristic frequencies. 

The subsequent analysis is concerned with methods of finding the roots 
of the infinite determinant. 

Although the present chapter does not include the effects of damping 
on rotor instability, Mr. George W. Brooks has prepared appendix B to 
indicate how effects of damping may be included in the analyses. 
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It is expected that designers will be able to obtain sufficiently 
accurate information by considering only the cases of infinite or zero 
support stiffness along one direction together with the cases of refer- 
ences 1 and 2. In order to avoid the necessity for extensive calcula- 
tions, a design chart is included giving the location of the principal 
self-excited- instability range for the case of infinite support stiffness 
in one direction. 


DERIVATION OF THE EQUATIONS OF MOTION 


The symbols used herein are defined in appendix A. 

The equations of motion are obtained from Lagrange's equations and 
from the expressions for kinetic and potential energy. Four degrees of 
freedom of the rotor system are considered: components of deflection of 

the rotor hub in the plane of rotation, and hinge deflections of the two 
rotor blades about their vertical hinges. All motions are thus assumed 
to occur in the plane of the rotor. The rotor is assumed to rotate at 
a constant angular velocity <x>. The analysis can be applied to rotors 
without hinges by assuming an effective spring stiffness and hinge posi- 
tion to represent the elastic deflection of the blade. 

The pertinent physical parameters are: 

a radial position of vertical hinge . 

b distance from vertical hinge to center of mass of blade 

mp mass of rotor blade 

m effective mass of rotor supports 

r radius of gyration of blade about center of mass 

K x , Ky spring constants of the rotor supports along the X- and 
Y-directions, respectively 

Kp spring constant of hinge self-centering spring 

Let the origin of the X,Y-coordinate system be placed at the undis- 
turbed position of the rotor hub. At time t equal to 0, the line 
through the blade hinges and rotor hub is assumed parallel to the X-axis. 
After a time interval t, let the rotor hub deflection be z and hinge 
deflections be (3p and respectively, where z is the complex posi- 

tion coordinate measured in a coordinate system rotating with the rotor. 
(See fig. 1.) Then the positions of the centers of mass of the two blades, 
as measured in fixed coordinates, will be, respectively, 
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The first term in equation (2) represents the kinetic energy of the rotor 
blades, including the energy due to rotation, and the second term is the 
contribution of the rotor hub. 

Upon expanding equation (l) into power series in Pi and £ 2 (only 
small deflections from equilibrium being considered) and substituting 
into equation (2), there is obtained 


T = | M(z + icuz) (z - ioa) + | m b -(p x - p 2 )(iu£bz ” ia;2bz + + ab2 ) 


(Pl - p 2 )(bil + cubz + abz - biz) + (h 2 + r 2 )(pi 2 + P2 2 ) * 

aba£(Pi 2 + p 2 2 )l ^ 


where only the terms that are quadratic in the variables have been 
retained, and M represents the total mass of the rotor system. 

The potential energy of the system V is given by 



^(z 2 e 2iajt + z 2 e“ 2iujt ) 


( 4 ) 


J 
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where K is the average stiffness and ZK is a measure of the dif- 
ference of the two principal stiffnesses; that is. 


K = 


Ky + K X 
2 


Z£ = 



As in references 1 and 2, simplifications in the analysis are intro- 
duced by replacing the hinge variables and P2 b y the following new 

variables : 


So - |(Pl + P2) 
01 - |(Sl - P2) 


In terms of these new variables the expressions (3) and (4) become, 
respectively, 


T 


^(z + io») (z 
2 


ioa) 


+ m b 



iu3z)(i0i - OJ 03 _) 


(z + icuz)(i9i + oc0^) + 



- | (eo 2 + e i 2 ) 


(5) 


and. 


V = — zz 


Kn ( 2 

+ -Il e 0 + 


ei 


2 ) _ ^f z 2 e 2ia3fc + z 2 e“ 2i ^t) 
4 \ > 


( 6 ) 
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By use of the Lagrangian form of the equations of motion 

_d /8t\ _ 8t + 8V = o 
dt\8x/ dx dx 

the following equations of motion for the rotor system are finally 
obtained: 

e 0 + (o£ai + A 2 )e 0 = o (7) 

(8) 
(9) 

o (10) 

where the notation 


dt 


(D + icu) 2 z + ip.(D + ico)^0i + — z - — ze"^-'- a ^ = 0 

M M 

(D - icu) 2 z - in(D - ia>) 2 6i + - z - — ze 2iajt = 0 

M M 

(D + iu>) 2 z - (D - ia)) 2 z + 2i ^1 + ^j^D 2 + Aia£ + A 2 j0;L 


has been used, and the following combinations of the original parameters 
have been introduced: 


H = 


2 % 

M 



A 2 


K P 
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Equations ( 7 ) to (10) do not include damping terms. The general 
equations for the damped case are derived in appendix B and it is shown 
that both the case for no damping and anisotropic supports, which is 
treated in the present paper, and the case for damping and isotropic 
supports, which was treated in reference 2, are readily obtained as 
special cases of the generalized equations by neglecting appropriate 
terms . 

Equation ( 7 ) can be solved independently of the others since it is 

an equation in 0 q alone. The equivalent equation appeared also in ref- 

erences 1 and 2, and its solution represents in-phase motion of the blades 
with no resultant reaction (except torsion) at the rotor hub. This motion 
will not be further considered herein. 

The problem is thus resolved into the solution of the three simul- 
taneous equations (8) to (10). It will be noted that the terms with 

periodic coefficients in equations (8) and ( 9 ) disappear if — =0, that 

M 

is, if Ky = Kx* Equations ( 8 ) and (9) are thus reduced to the problem 
treated in reference 2. 


FORM OF SOLUTION OF EQUATIONS OF MOTION 


The equations of motion (eqs. (8) to (10)) are similar in mathematical 
properties to Mathieu's equation, which occurs in the analysis of vibrating 
systems of one degree of freedom with variable elasticity. (See ref. 4.) 

A generalized form of Mathieu's equation was solved analytically by Hill. 
(See ref. PP« 413-417*) An extension of Hill's method has been applied 
in reference 3 to a problem involving two degrees of freedom, and a fur- 
ther development of the method of reference 3 is followed in the present 
paper . 


Equations (8) to (10) constitute a system of linear differential 
equations with periodic coefficients. Three second-order equations possess 
six linearly independent solutions that, according to the Floquet theory 
(ref. 5 , p. 412), are of the form of an exponential factor times a periodic 
function of time. Particular solutions are of the form 


Z = e^Ptt) + e" ia) a t Q(t) 
z = ei^ Q ( t ) + e -i<%tp( t ) l 

0]_ = e _: *- ( %^R(t) + e - ^%^R(t) 


( 11 ) 
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where c% is known as the characteristic exponent, and P(t), Q(t), 

and R(t) are periodic functions of period n/ca. 

Since P(t), Q(t), and R(t) are periodic functions of t, they 

can be represented by complex Fourier series, and equations (ll) become 


= Jr A je^ 2la>Hj;) a) it + 


£ = Y_ + ^T Aie -(2lcDK%)it 


( 12 ) 


e l = J Cl e( 2lu>H %) it + c^ 2 * 0 * 0 *^ 


where A^, B^, and are complex constants. 

Equations (ll) and (12) show that, when the rotor system is stable 
and cua is real, the motion not only is not simple harmonic as was the 
case in references 1 and 2, but, in general, is not even periodic. The 
motion can be said to consist of a fundamental frequency c% plus 

"harmonics" of frequency on, + 22 cd where 2 is any integer. From equa- 
tions (12) it is seen that the value of c% is not uniquely determinate, 
since + 22o> also satisfies equations (12). (The imaginary part of 
cca is definite, however.) It can be shown furthermore that, corresponding 
to each value of oca, -oa is also a solution. Only those three values 

of oa for which the real parts lie between 0 and <n need therefore 

be considered. These values will be referred to hereinafter as the three 

"principal" values of oa- These values of differing from the prin- 

cipal value by 22cd, or having opposite sign, will be referred to as 
"harmonics" of the corresponding principal value. 

Since z has been defined as a position coordinate in a rotating 
frame of reference, the values of oa can be interpreted as the natural 

frequencies of the rotor system in rotating coordinates. 
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SOLUTION OF EQUATIONS OF MOTION 


Determinants! Equation 


If the formal solution (eqs. (12)) is combined with the equations 
of motion (eqs. (8) to (10)), and the coefficients of each exponential 
time factor is separately equated to zero, an infinite set of homogeneous 
equations is obtained. These equations can be separated into two inde- 
pendent sets. Each equation of one set is the conjugate of an equation 
of the other set, and only one set need be considered. Thus 


< | - [ia + (21 + l)w] 2 l A l - Bj +1 - (lijofe + (21 + 1)^1% = 0 

1 ' ( 13 ) 


v r -2 

c% + (21 - l) o 

M L -1 


Bz - ^ A 
L M 


2 

^1-1 + + (21 - l)o>] C j = 


(14) 


-[fife + (21 + 1 )od] Aj + |ua + (21 - 1)cdJ + 2i /l + 


-( c % + 21cd) + 


Aqo? 



C i 


= 0 


( 15 ) 


where 1 takes on all’ integral values from -« to <». 

In order that the values of Aj, B j, and Cj not equal zero, the 
determinant of the coefficients of Aj, Bj, and must be zero. This 

determinantal equation is 
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a.4,-4 a -U,-3 000 

0 

0 

0 

0 

a -5,.4 a -3,-3 a -3,-2 0 0 

0 

0 

0 

0 

0 a -2,-3 a - 2 ,-2 a -2,-l 0 

0 

0 

0 

0 

0 0 a -l,-2 a -l,-l a -l,0 

0 

0 

0 

0 

0 0 0 ®0,-l ®0,0 

*0,1 

0 

0 

0 

0 0 0 0 a 1 ^Q 

a l,l 

a l,2 

0 

0 

0 0 0 0 0 

* 2,1 

* 2,2 

*2,3 

0 

o 

o 

o 

o 

o 

0 

a 3,2 

a 3,3 

a 3A 

0 0 0 0 0 

0 

0 

%3 

a M 


= 4 <%) = o 


(16) 


where 


i -4,-4 


= -1 + 


k/m 

(<x& - 3oo) 2 


*o,o - - 1 + 


A]_cu4 + A 2 




1 -4,-3 = a 5 




+ CD) 


/ 0 ^ - 3 cd 
^-3,-4 “ l OCfeL - 2CD y 


a l,0 - a 3 
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-3,-3 


a-ju£ + A2 
(<x& - 2a>) 2 


1,1 ’ - 1 + 7^7 

( cife . + Co ) 


co& - co 
<% - 2“ 


zk/m 

11,2 (<% + co) 2 


l -2,-3 - a 3 


zk/m 

l 2,l - ; “i 

(c% + co) 


T . **■/ * * 

a_p -2 = + o 

(<% - m ) 2 


= -1 + 


(oca. + co)‘ 


zsc/m 

l -2, -1 ■ p 

(coa - co) 


a 2,3 = A 3 


l ' 1,_2 (<% - co) 2 


c % + CO 
afe, + 2o) 


a -l,-l = - 1 + 


(c% - co)‘ 


a 3,3 = _1 + 


Aj(i£ + Aq 

r 

(c% + 2m) C 


l -l,0 = A 3 


/(% + 3tO 

^ ,hr \c% + 2 m, 


a O,-l 


c% " co 


a 4,3 " A 3 

, k/m 

a 4 4 ~ -1 + 2 

(a* + 3co) 2 
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The determinant has been somewhat simplified by multiplying and dividing 
the rows and columns by various quantities, and the parameter has 


been substituted for its equivalent 


n/2 

1 + (r/b) 2 


Let this infinite determinant be A((%) . The problem consists in 
solving the equation A(<%) = 0 for its roots c%. These roots will be 
infinite in number, consisting of the three principal values of plus 


all their harmonics. The values of 




as a function of 


' 0 u/\l k/m. 


are seen to depend only on the values of the three nondimens ional par am - 


A ^ 

eters An , — and At: and the stiffness -ratio parameter AC/k. 

k/m j 


A determinant of infinite order has meaning only insofar as it is 
defined as the limit of a determinant of finite order. Define ^(o^) 

as the determinant of order 6n - 3 formed from a square array of A(c%) 

An oo2 + A o 

centered on the term -1 + — This term, which originally was 

<% 2 

associated with Co in equation ( 15 ), will be referred to hereinafter 
as the "origin" of the infinite determinant A(c%). The choice of this 
term as center of A((%) is purely arbitrary, and it was selected solely 
for reasons of symmetry. 

Thus 


A(<%) = lim Ai(<ha) 
n — >oo 


(17) 


The limiting values of the roots of the equation 


4i(<ha) = 0 


as n becomes infinite will be the values of the roots of the equation 


A(<jQa) = 0 



IJ 
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The method of calculating the roots of A(oa) = 0, by successively 
calculating the roots of Ai(ufa.) = 0 for larger values of n, is entirely 

too tedious. Instead, the method reference 3 will be followed. This 
method involves the calculation of the value of a(c%) for several spec- 
ified values of ofe. The roots of A(a&) = 0 can then be obtained from 
a trigonometric equation involving the roots and the calculated values 
of A(t%) . 


Auxiliary Determinants and Recurrence Relations 
for Calculating A(o&) 


Before the trigonometric equation is derived, it is convenient to 
have a systematic numerical procedure for determining the value of A(c%) . 

As . n becomes infinite, the terms of extend to infinity both 

above and below the origin. By expanding in terms of the ele- ‘ 

ments of the column containing the origin, it can be expressed in terms 
of auxiliary determinants that extend to infinity in only one direction. 
Recurrence relations can then be obtained that give the value of these 
auxiliary determinants . 

The auxiliary determinants are minors of ^(ofe) and are defined 
as follows : 


Cn(°h.) 


determinant of order 3^-2 consisting of the terms below 
and to the right of the origin; that is, determinant having 


first row and column beginning with term 


i + — s/a — 

(<% + (jd ) 2 


Dn(<%) 

En(cua) 

Ln(tba) 

Mn(<ha) 

Nn(ote) 


determinant 

of 

order 

3n - 

3 

last row 

and 

column 



determinant 

of 

order 

3n - 

4 

last row 

and 

column 



determinant 

of 

order 

3n - 

3 

first row 

and column 


determinant 

of 

order 

3n - 

4 

last row 

and 

column 



determinant 

of 

order 

3n - 

5 

last row 

and 

column 




formed from 

C n M 

by 

omitting 

formed from 

Dn (<%) 

by 

omitting 

formed from 

Cn(<Ha) 

by 

omitting 

formed from 

Ln(tba) 

by 

omitting 

formed from 

M n (<ha) 

by 

omitting 
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The following three determinants will also be needed: 


Gn(<%) 

determinant of order 3n 
first row and column 

- 4 

formed from 

Ln(tth) 

by omitting 

Hn(<%) 

determinant of order 3n 
last row and column 

- 5 

formed from 

Gn(“a.) 

by omitting 

In^) 

determinant of order $n 
last row and column 

- 6 

formed from 

H n (<%) 

by omitting 


Determinants similar to the foregoing can be formed in the same 
manner from the upper half of £n( c%) • Denote these determinants by the 

subscript -n instead of n. It is seen, however, that their values can 
be obtained from the values of the determinants already defined from the 
lower half of by merely replacing with -oca (for example, 

C-n(<ra) = Cn(-oa)). 

• Expanding Ai(<%) in terms of the elements of the column containing 
the origin gives 


Ano^ + Ap\ 

£*(<%) = 1-1 + — -|Cn(c%)C n (-uja) 


(tier 


a 3 


j L n (-o%)C n (oa) 


+ 


f \2 

— -)Ln(<%)0n(-aa) 

V / 


(18) 


The auxiliary determinants C n (c%), Dn(oa,), a 11 *! En(<ha) satisfy 
the following recurrence relations (obtained by expanding each in terms 
of the elements of its last row) : 
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Cn(t%) = 


-1 + 


K/M 

[(% + (2n - 1 )gd] 2 




■ t% + (2n - 1) 
+ (2n - 2) 


00 


00 


En(ufe) 


Dn(aa) 


= 1-1 + AjO^ + A2 1 


®n( aa) — 


-1 + 


[oa + (2n - 2)cp] 2 

K/M 

[(Ob. + ( 2 n - 3 )co| 2 1 


En(<%) - A3 


-n-l(<%) 


<% + (2n - 3)ao 
j% + (2n - 2)cd 

[oca + (2n - 3 )<jd ]4 


Cn-l(ofe) > (19) 


Dn_i(o)a) 


The determinants L^oa), M n (aca)> and N n (ua) and also the system 
GnCafe.)* Hn( oa ) > and In(aa) satisfy the same recurrence relations as 
C n (<hab Dn(aa)^ and E n (oa) , respectively. 

The values of these nine determinants can be found from the recur- 
rence relations (eqs. (19)) and the following initial values, obtained 
directly from equation (l6) : 


c l(<%) = -1 + 


E 2 (c%) = 


N 2 (<%) = 


M2(<%) = 


H2(aa) = 


G 2 (ah.) = 


k/m 


(<% + co) 2 



1 1 K/M 

2 

(a k/m) 2 


(<% + a>) 2 


(<% + co) 1 ^ 


1 , . k/m 




(<% + to) 2 




r~ — 

! K/M 

i _ 

A^ocf- + A 2 


(<% + (o) 2 

JL — 

((% + 2cd) 2 


A-i oc? + Ap 

1 i x . . c 




(aa + 2 cd) 2 




r + A 2 

IL 

k/m 


- A31 


^c% + co N 

+ 2a V 


(<% + 2ai) 2 


(<% + 5 a)) 2 


A 3 


+ 3 <x> 
{<Oa. + 2 <Oj 


,2 


( 20 ) 
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By use of the initial conditions (eqs. (20)), the recurrence rela- 
tions (eqs. (19)), and equation (l8), the value of Aa(a&) can be cal- 
culated. The value of A(o&) will then be. the limiting value of £n(ofe) 
as n becomes infinite. 


The Behavior of ^(ota) for Large Values of n 

So far it has been tacitly assumed that the determinant A(<xa), as 
defined in equations (16) and (17), is convergent, and further, that it 
remains a function of aa in the limit as n becomes infinite; that 
is, it is not identically equal to zero, independent of the value of c%. 
It will now be shown that the function Ai(aa) does become zero in the 
limit j independent of <%, but that when Z^(aa) is divided by an appro- 
priate function of n, a new function F n (ua) will be obtained which will 
be convergent and remain an unambiguous function of in the limit. 

The derivation of the appropriate function by which to divide ^(oa) 
evidently depends upon the behavior of £a(aa) as n becomes very large. 
As n becomes infinite, the recurrence relations (eqs. (19)) become 


On = -Dn - A^En (21a) 

Dji = -En - AjC n _i (21b) 

En = -Cn-1 (21c) 

with identical equations for Ln, Mn, and N n and for Gn, H n , and In- 

Equations (21) are readily solvable since they constitute a system of dif- 
ference equations with constant coefficients. They are satisfied by solu- 
tions of the form 


C n = C 0 k n (22a) 

whdre Co is some arbitrary constant and k is a constant to be deter- 
mined. From equations (21c) and (21b), respectively, 

E n = -CQk 11 - 1 


(22b) 
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and. 


Dn = Cok 11 " 1 - AjCgk 11 - 1 


= Cq(1 - 


(22c) 


Combining equations (22a), (22b), and (22c) with equation (21a) and 
dividing through by CQk n- l gives 


k = 2A 3 - 1 (23) 


Thus, by use of equation (l8), it is seen that for large values of n, 
Ai(c%) varies as k^ 11 . Since by definition Aj must have a value 

between 0 and l/2, k^ must lie between 0 and 1. Thus, in the limit, 


lim Ai(<x&.) = 0 


independent of c%. 

Consider the function 


F n (<ba) 


k 2n 


The equation F n (o^,) = 0 will obviously have the same roots for 
as does £n(<%.) = 0. The function F n (<%) has the advantage, however, 

as is seen from the preceding discussion, of remaining an unambiguous 
function of in the limit as n becomes infinite. Define this 

limit as 


F(ofe) = lim F n (<ra) 


n- 


= lim 


Ai(oh) 

k 2n 


n- 


(24) 
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The primary problem can now be redefined as the problem of determining 

the roots, infinite in number and consisting of <%_, , and 

12 3 

and all their harmonics, of the equation 


F(afe) = 0 


(25) 


Evaluation of Roots of Equation (25) 

The following trigonometric expression for F(c«a) will now be 
derived: 


F(<%) 


lim 
n — 


Ai(<%) 

k 211 



(26) 


The function F((%) 
to be periodic of period 


is seen from equation (24) and equation (l6) 
2 cd, to have roots * 2sooj, *( a te > 2 * 


and t^aco.^ - 2soo) 


where 


s is any integer, to have second-order poles 


at ofe = ±2 s<d, and to have fourth-order poles at = ±(s + l)a). 
Liouville's function theorem states that a function of a complex variable 
( in this case, ace.) that is analytic everywhere in the complex plane, 
including the region at infinity, must be a constant. It will be shown 
that F(cife) is finite at infinity (except if <% proceeds to infinity 
along the real axis). If the poles along the real axis could be elim- 
inated by forming a suitable function of F(o)a)> without at the same time 
introducing new poles, then that function, by Liouville's theorem, must be 
a constant. 


Such a function of F(ac&), which is analytic everywhere in the com- 
plex plane, is 


J(c%) = 




\2co ) 


5 

IT 

0=1 


lr.2^ - sin2| 

\2coj 


sin* 


, 2co 


(27) 
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where cua-^> <xq. 2 , and are the three principal values of c%. The 

function J(c%) is therefore a constant. The value of j(cqa.) found by 
making c% approach infinity along the imaginary axis is 

j(afe) = P(») 

where F(<») is the value of F (on.), as o^, becomes infinite. 

The value of F(°°) can be found by letting in ^(0%) 

and then letting n — >00. From the form of Z^(~) it follows that 

4i(°°) = NgnC 00 ) = -C2n(°°) 

The recurrence relations defining C n (°°) are the same as equations (21). 
The expressions Dn and En may be eliminated from equations (21 ) , and, 
therefore, C n may be given as follows: 

C n = -(1 - 2A 3 )C n _i 

The initial conditions (eqs. (20)) reduce to = -1, from which it 
follows immediately that 


Z^(co) = -(1 - 2A3) 211 " 1 


= k 211 " 1 


from equation (25)- Therefore 


F(oo) 


= lim 


4i(”) 

k 2n 


lim 
n — >oo 


k 211 - 1 

k 2n 



(28) 


Thus equations (27) and (28) lead to the evaluation of equation (26). 




After equation 
oqa 2 , and 


(26) bas been obtained, the problem of determining 
may be considered theoretically complete inasmuch 
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as equation (26) is really an identity in o^. Suppose that o& is 
assigned any specific value in equation (26) and F(c%) is computed to 

a certain degree of accuracy. If these computations are made for two 
more values of all different, equation (26) will have yielded three 

equations in the three unknowns u^, <rn 2 , and u^y These equations 

can then be solved for the principal values of Any degree of accuracy 

may be achieved by carrying out the computations for F(c%) to a suffi- 
ently large value of n. 

The foregoing procedure can be systematized by rewriting equation (26) 
as 





(29) 


A convenient choice for the three arbitrary values of is o^. = 0, 

(jd, and <d/ 2. The explicit definitions of K(0), K(l), and K(l/2) then 

become 



The equations for evaluating K(0), K(l), and K(l/2) are 


3 J 
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^<^> sln2 (lir) oos4 (5?) 

kF(c%) sin^-^j cos^ 


K(l/2) = lira 

c% — W 2 


Carrying out the limiting operations indicated in equations ( 31) and 
using the auxiliary determinants C n (c%), Dn(o&), and so forth, give 





K(0) = lira 


K(l) = lira 
>oj 


K(0) = lim 
11 


j^A;ja£ + C n ^(o) - 2A^a£l^(0) C n (0) 


-*ko 2 (l - 2A 5 ) 


2n-l 


K(l). = 


lim 
n — >oo 


« 4 £ ^ G n 2 (-a>) 
M M 

_16(1 - 2A 3 ) 2n ' 2 _ 


(32) 


K(l/2) = 


4(a]o£ + Ag) 

Cn (; 

Ht) - 

Ln 

£2>W“) 

^2 J \2) 

* *#“(?)] 

L J 


-8(1 - 2A5) 211 ' 1 

L 




where the quantities in brackets are conveniently represented by K(0) n , 
K(l) n , and K(l/2) n , respectively. The quantities K( ) n are used in 
numerical computations as approximations to the functions K( ) . 

The formulas (32) for K(0), K(l), and K(l/2) converge slowly 

with increasing n. The convergence can be speeded up greatly by making 
use of the concept of convergence factors used in reference 3* A con- 
vergence factor for K( ) n is a function of n approaching the limit 1 
as n becomes infinite, which, when multiplied by K( ) n , gives an expres 
sion which converges' rapidly with increasing values of n to the values 
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of K( ) . The details of the derivation of an appropriate convergence 
factor for K(0) n will be found in appendix C. Convergence factors 
for K(l) n and K(l/2) n are derived in a similar fashion. The resultant 

expressions are 


K(0) = 


lim 
n — ^oo 


K(0) n sin 2 ^ 

* to - 


4 j=l 


4J 


,2 


K(l) = 


lim 
n — 


K(l) n cos 2 


Q 


Q 


n-1 / 

TTf 1 - - 

j=l \ (2j - ifj 


(33) 


K(l/2) = 


lim 


K(l/2) n cos 


2n-l 

"TT 




4Q \ 
(2j - if) 


y 


where 


Q = 


2 ^(1 - Aj) + Ayx£ + A2 + 2A 50 ^ 


o^(l - 2A3) 


For a given value of n, the quantities in brackets are found to be 
better approximations to the respective values of K( ) than K( ) n 

alone . 

The method of obtaining the values of may be summarized as 

follows. By vise of the initial conditions (eqs. (20)) and the recurrence 
relations (eqs. (19)), the values of the determinants C n (0), Ln(0), 

Gn( -cd) , C n (u72), C n (V 2 )> Ln(^/2),and Ln(-o>/2) can be computed for 
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increasing values of n. With the substitution of these values into 
equations (33) and with the use of equations (32), approximate values 
of K(0), K(l), and K(l/2) can be computed. The process appearsto 

be rapidly convergent with n, especially for large values of ay|rK/M. 

The values of cjc^, anc ^" can ^ en found, from equations ( 31) , 

the definitions of K(0), K(l), and K(l/2). 


Conditions for Stability 


From equation (13) the condition for stability of the system is seen 
to be that all three values of must be real numbers. If any one of 

them is complex or pure imaginary, then one of the terms in the solution 
(eqs . . (ll) ) will increase indefinitely with the time, the motion therefore 


being unstable. 


This condition implies that the expressions 


sin^l 


rjtc%2 

, 2cd 


sin 2 £ I and sin 2 all are real positive numbers less than or 

\ 2cd I ’ \ 20) / 

equal to 1. The conditions for stability can be expressed directly in 
terms of K(o), K(l), and K(l/2) by means of their definitions 

(eqs. (30)). The three equations (30) are formally equivalent to a 
single cubic equation 


x3 + bx2 + cx + d = 0 


the roots x-j_. 


x 2 , and x^ of which are 




and 


so forth, and the coefficients b, c, and d of which are functions 
of K(0), K(l), and K(l/2) where 


2b = 4K(0) + 4K(l) - 8K^j - 3 
2c = -6 k(0) - 2K(1) + 8k/|) + 1 


d = K(0) 
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After some manipulations involving the Descartes rule of signs, the nec- 
essary and sufficient conditions for stability are found to be 


0 ^ -K(0) % 1 
0 1 K(l) ^ 1 
-1 < 8K(l/2) < 1 


A 


> (34) 


A = l8bcd - 4b3d + b 2 c 2 - 4c3 - 27d3 > 0 

y 

The quantity A is the discriminant of the cubic equation. 


SPECIAL CASES OF GENERAL THEORY 


Three special cases of the general theory are of interest. These 
cases are the cases for which one of the principal stiffnesses Ky is 

respectively zero, equal, or infinite in magnitude in comparison with 
the second principal stiffness K x . 


Case of Ky = K x 


The case of Ky = K x has been treated in reference 2. If Ky = K x , 

equations (8) to (10) reduce to the equations of reference 2. In this 
special case, the motion of the rotor system becomes simple harmonic, 
since all the coefficients A^, B^, and C^ in equations (12) are 

identically zero except A q, Bq, and Cq. 

Case of Ky = 0 

The special limiting case of Ky = 0 is of interest in the case of 

a pylon of which the stiffness is negligible along one principal direction 
with interest centered on the frequencies involving the other principal 
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stiffness. In the case of Ky = 0, the function K(l) as given by equa- 
tions (32) becomes identically zero. This result is also evident from 
the original definition of K(l) as given in equation (28), because one 
of the values of o^, say a is of necessity equal to +o>. (it will 

be recalled that ofe is the frequency as measured in rotating coordinates. 
In fixed coordinates it would be zero.) 


It is possible to give much simpler stability criterions for this 
case because there are only two K-functions, K(O) and K(l/2), and 
two values of d%, and (%£, to be determined. The new K-functions 

may be defined as follows: 


Kl =. 



k 2 = 




-K(0) 


> 

1 + 2K(0) - 8k(|J 
2 


(35) 


In terms of Kp and K2, the criterions for stability become 


0 < Kp ^ 1 

(36a) 

0 ^ K 2 $ 1 

(36b) 

+ ^2 ^ 1 

(36c) 


Given the values of Kp and K2, the values of o^p and can 

be determined from equations (35)* A graph of the relation in equa- 
tions (35) is given in figure 2 by means of which the real values of 

and can be read off directly once Kp and K2 are known. 

o 

A graph showing the variation of Kp and K2 with , for 

K x /M 

the typical parameters Ap = 0.1, A2 =0, A3 = 0.1, and Ky = 0, is 
shown in figure 3- By use of figure 2 the values of and c% 2 


can 
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be obtained. These values are shown in figure 4 plotted against oy^K x ^M. 
Calculations are carried down to 03 = 0.5- The general behavior below 

VM^ 

this speed is discussed in the section entitled "General Behavior of Rotor 
System as a Function of Rotor Speed." 

Case of Ky = oo 


The formulas for the limiting case of Ky = °° cannot be obtained 

conveniently from the general theory. Instead of carrying out the limiting 
process, it appears preferable to begin by treating the problem as one of 
only three degrees of freedom (two hinge-deflection coordinates and one 
hub-position coordinate x) and by developing the theory along lines 
similar to those used for the general treatment. In this way a system 
of two simultaneous equations with periodic coefficients is obtained, 
with the variables 0q and x. These equations are solved in a manner 

similar to that for the general case, the treatment being simpler, how- 
ever, since the solution has only two principal values of <x&. 


The details of the solution of the equations of motion, together 
with the final formulas for the K-f unctions, including convergence fac- 
tors, are given in appendix D. It is found that the same Kq and K2 

occur as for Ky = 0. The criterions for stability are exactly the same 
as those for Ky = 0, the conditions of equations (36). Figure 2 can 
also be used to determine the values of from the values of Kq 

and K2 • 


A graph giving the variation of Kq and K2 


parameters Aq = 0.1, A2 = 0, A^ 
ure 5- In figure 6 the values of 
plotted against o)A/K x /M. 


= 0.1, and 

<%qtec/M 


K-u 


and 


with 


K x /M 


for the 


= 00 is shown in fig- 




^1 


M are shown 


DISCUSSION OF RESULTS 


Types of Instability 

Instability may occur as a result of the violation of any one of 
the stability criterions of equations (34). Violation of each condition 
is associated with a different type of instability, which would show up 


NACA TN 3844 


95 


differently in the motion of the rotor system. Experience with computa- 
tions indicates, however, that the criterions of most practical importance 
for helicopters are 


A > 0 
-K(0) ^ 0 
K(l) ^ 0 


Similarly, the important criterions in the limiting cases of Ky = 0 
and Ky = «> are 


/Kl + \/K2 .= 1 


K x ^ 0 
K 2 ^ 0 


If the condition A = 0 or + \/K2 = l) is violated, the other 

conditions being satisfied, then and will be complex con- 

jugates, and the rotor system will execute self -excited vibrations at 
frequencies, in general, incommensurate with the rotor speed. (Higher 
har monics will also be present.) This type of instability will be 
referred to hereinafter as a "self-excited vibration." 

If the stability condition -K(0) ^ 0 (or Kq ^ 0) alone is violated, 
then one of the values of <x& will be a pure imaginary number. Physically, 
the rotor system will execute self-excited vibrations having a basic fre- 
quency, as seen in rotating coordinates, of zero. This behavior is sim- 
ilar to the ordinary critical-speed behavior of a shaft. Frequencies at 
higher harmonics 2nu> will also be present. This type of instability 
will be referred to as a "self-excited whirling." 

The third stability condition K(l) =0 cannot be violated since 
K(l) as given by equation (32) cannot be negative. However, K(l) can 
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be exactly equal to zero. (A similar statement applies to K 2 O At such . 
a point, where the rotor system is on the border line between stability 
and instability, one of the values of will be equal to ±<n. In 

fixed coordinates this result means that the rotor system will have a 
natural frequency equal to zero. The rotor system will, therefore, be 
in resonance with a steady force - a force of zero frequency. The ampli- 
tude of the zero-frequency term for the hub motion in such a situation 
can be shown to be zero, but the blades will oscillate. Also, higher 
harmonic terms, notably the term of frequency 2cd (in nonrotating' 
coordinates), will show up in the hub motion. This type of vibration, 
which is a resonance phenomenon and not a self-excited vibration, will 
be c all ed a "steady- force resonance" vibration. 

Each of the vibrations described - self-excited vibrations, self- 
excited whirling, and a steady-force resonance vibration - appeared in 
the discussion of the two-blade rotor on equal supports (ref. 2); how- 
ever, there the motions were simple harmonic, no higher harmonics being 
present. - 

General Behavior of Rotor System as a 
Function of Rotor Speed 

The approximate location of the instability regions can easily be 
found by examining the limiting case of Aj = 0, that is, the case of 

zero coupling between the blade and hub motions. For simplicity, the 
discussion is also restricted to the case of free hinges (A 2 = 0) and 

Ky = 00 . The Kq and K 2 functions become 

cos 2 (— 

\2cjd 

K 2 = 

Eliminating the rotor speed cn from equations (57) gives 

K 1 + — -= 1 . 

cos2 (ffl) 
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■which considered as an equation in the variables Kp and K 2 represents 

a straight-line segment (one of the lines in fig. 2) terminated by the 
Kp and K 2 axes. The segment can be shown to be tangent to the curve 

^ + y*i= 1 . as a) decreases, the representative point moves up and 

down the line segment, performing an infinite number of such oscilla- 
tions as co approaches zero. Whenever Kp = 0, the point is at a self- 
excited- whir ling speed. The corresponding speed is 


\[Kx/M 

2s + 


where s represents any positive integer. Thus a self -excited whirling 
will occur when the rotor speed is approximately equal to 1, l/3, l/5> 

l/7, and so forth of the natural frequency of the hub Similarly 

it can be shown that there will be a steady-force resonance vibration when- 
ever the rotor speed is approximately equal to l/2, l/4, l/6, and so forth 

of the hub natural frequency \Jk x /M. Finally self-excited vibrations will 
occur at rotor speeds approximately equal to 


CD = 


VKx/M 

•2s + 1 - 


Figure 7 shows the general pattern of response frequency plotted 
against rotor speed for a small value of the mass-ratio parameter Aj. 

Along the horizontal parts of the curves, blade motion predominates over 
pylon motion. Pylon motion predominates along the slanting parts. 

Although the foregoing discussion was developed for the case of 
Ky = co, it is believed to apply equally well to the case of Ky = 0 

and also to the general case of K y =/ K x if the ro ' tor hub is considered 
to have two natural frequencies \jK x jM and \jKyjlA, each frequency having 

associated with it an infinite set of instability ranges located at approxi 
mately the speeds given. 
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Comparison of Results for Different 


Values of 



Figures 4 and 6 give the principal values of If] 

d^Kx/m 


against the rotor speed to/ 


for 


Ky = 0 


M plotted 


and K v = «>, respectively, 


both calculated for the same set of parameters Aq, A2, 
calculations have been carried down to - = — = 0.4. 




and Aj . The 
The similarity 


between the two curves is striking. So far as the calculations have been 
carried, each system shows the presence of one self-excited-vibration 
instability range, one self -excited-whirling instability range, and one 
steady-force resonance speed A. If the calculations were carried to 
lower values of u>, further instability ranges and steady-force resonance 
speeds would appear. 


For comparison, the response frequencies of a two-blade rotor on 
equal supports (Ky = K x ) for the same set of parameters is shown in fig- 
ure 8 . The frequencies were calculated from the formula in reference 2. 

Down to 03 - = 0.5, this chart is very similar to figures 4 and 6. In 

l/Kx/M 

addition it shows one range of rotational speed at which self -excited- 
vibration instability occurs, one range of rotational speed at which 
self -excited- whirling instability occurs, and one range of rotational 
speed at which a steady-force resonance speed occurs. Figure 8 differs 
principally from the figures for Ky ^ K x in that it shows no further 

instability ranges at low values of od. 

In references 1 and 2 charts are presented giving the location of 
the self-excited-vibration instability range for various values of the 
parameters Aq, A 2 > and A^. A similar chart for the case of a two- 

blade rotor with Ky = °° is given in figure 9 . In using the chart, a 

straight line is drawn representing the variation of — — with the 

K x/ M 

function Ala ^ — - A2 . The intersections of this line with the appro- 

K x /M 

priate Aj curves give the beginning and end points of the instability 

range. The dashed line in figure 9 indicates the stability ranges for 
the parameters of figure 6 . 

Some observations concerning the relative location and extent of 
the various instability ranges in figures 4, 6, and 8 appear to be appli- 
cable to a wide range of values of the parameters Aq, A 2 , and. Aj. Thus 


NACA TN 3844 


99 


the self -excited-vibration instability range in the case of Ky = K x 

(fig. 8) is wider (and hence the vibration probably more severe) than 
the corresponding ranges in the cases of Ky = 0 and Ky = <». (See 

figs. 4 and 6.) Also, this instability range occurs at lower rotor 
speeds in the case of Ky = 00 than it does in the cases of Ky = K x 

and Ky = 0. The self -excited-whir ling instability range is considerably 
narrower for Ky = °o than it is for the Ky = K x case, and it is still 
narrower in the Ky = 0 case. 

In the general case of Ky 4 K x the location and extent of the 

instability ranges can be found fairly accurately by considering the 
problem ‘as the superposition of two problems, one of finding the signifi- 
cant rotor speeds referred to ^K X ^M as reference frequency with Ky 
assumed infinite and the other of finding the significant rotor speeds 
referred to yKy/M as reference frequency with K x assumed zero. With 

the foregoing discussion as a guide, sufficiently accurate design informa- 
tion can be obtained without extensive calculations for each value of 
KyjK x encountered in practice. 


Effect of Damping 

Although the effect of damping has not been examined mathematically, 
because complications would be introduced in the analysis, several infer- 
ences from the damping investigations in references 1 and 2 can probably 
be safely applied to the rotor-system studies in the present report. The 
numerous instability ranges occurring at low rotor speeds, which are very 
narrow and represent a mild type of instability, are probably completely 
eliminated by the presence of a slight amount of damping in the rotor 
system. The primary self -excited-vibration instability range can probably 
be narrowed and eliminated by introducing sufficient damping into both 
the rotor supports and the blade hinges. 


APPLICATION TO DUAL ROTORS 


It is easily shown that the analysis for the case of Ky = <» applies 

also to the case of a counterrotating rotor system consisting of two equal 
two-blade rotors revolving at equal speeds and acting equally upon the 
same flexible member. The rotors may be on the same shaft or on different 
shafts so long as the nonrotating flexible member is the same for both 
rotors. The supports, moreover, may have unequal stiffness in the X- 
and in the Y-directions, provided that the undeflected blade positions 
make equal angles with a principal stiffness axis. 
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The proof consists in showing that the energy expressions for the 
dual-rotating system can be separated into two independent sets of terms, 
each of which is of the same form as for a single rotor with K y = «. 

The resulting equations of motion will thus also be the same . 

The separation is accomplished by introducing new variables 


£l “ 2 ( 01 P° S ” 01 neg) 

= 2( 9 °P os " 9 ° ne s) 

and 

^1= !( 0 lpos + 0 lneg) 
^0 = g^Opos + 9 °neg j 


where the subscripts pos and neg refer to the 9 s defined for the 
rotor turning in the positive direction and for the rotor turning in the 
opposite direction, respectively. The energy expressions become 


T = 


i Mx 2 + 2m b -2x^ii sin ^ cos a3fc ) + 

(x ♦ ♦ tf) - i ^ + ‘o 2 )] * 

2y^r|i cos 0* - sin + 


| My 2 + 2m b 


2K, 


V = 


-^l 2 + ^o 2 + ^l 2 + ^o 2 ) + | Kxx2 + | 


(58) 


where M is the total mass of the system (M = m + 4)%) • 
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From equations (38) it is seen that | is coupled only with x 
and T) is coupled only with y. Equations (38) yield equations of motion 
of the same form as equations (D2) and (D3)« 


The stability properties for the dual-rotating case are thus exactly 
the same as in the case of Ky = °° for the single-rotating two-blade 

rotor. In particular, figure 9 can be used to find the location of the 
primary self-excited- vibration instability range. The value of A3 for 

2m b 

the dual-rotating rotor is defined as A3 = rather 

(- + H) (1 ♦ $) 

m b 

than A3 = as for the single -rotating rotor. All other 

( m + * $) 

parameters are the same for both cases. 


The quantity sin cut can be interpreted physically as the 

x-component of the displacement of the center of gravity of the blades 
due to hinge deflections. The quantity cos ojt is the corresponding 

y-component. The separation of the variables means, physically, that the 
motion of the system can be separated into two independent modes, each 
of which involves linear motion of the supports along one of the prin- 
cipal stiffness axes. 

Similarly, the stability of counterrotating rotor systems of six or 
more equal blades can be determined from the results of reference 1 with 

Ky = CO. 


CONCLUSIONS 


The following conclusions are indicated by the results of an investi- 
gation of the problem of vibration of a two-blade helicopter rotor on 
supports that have different stiffnesses along the two principal stiffness 
axes : 


1. Many speed ranges are found in which self -excited oscillations 
can occur. These oscillations are of two types - self-excited vibration 
and self -excited whirling. There are also many speeds at which steady- 
force-resonance vibration may occur. 

2. A stability chart which shows the self -excited vibration, self- 
excited whirling, and steady-force resonance speeds of the highest rotor 
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speed for each support natural frequency for a two-blade rotor on aniso- 
tropic supports is s imi lar in appearance to a stability chart for a two- 
blade rotor on isotropic supports. However, for the same rotor parameters, 
the instability regions are changed somewhat in position and extent. 

3* Mild self-excited-whirling speed ranges exist at rotor speeds 
approximately 1 / 3 , l/5> l/7> and so forth of each support natural fre- 
quency. Steady-force resonance speeds exist at approximately l/2, l/k, 
l/6, and so forth of each support frequency. Self-excited vibrations 
also occur at certain low rotor speeds. All these mild instability ranges 
are probably eliminated by the presence of moderate amounts of damping 
in the system. 

it-. A familiarity with typical results of limiting cases of the 
support- spring constants Ky = », Ky = K x , and Ky = 0 should enable 

a designer to avoid extensive calculations of cases of unequal support 
stiffness. In the general case of unequal support stiffness, the loca- 
tion and extent of the instability ranges can be found fairly accurately 
by considering the problem as the superposition of two problems, one of 

finding significant rotor speeds referred to one support frequency \/k x /m 
as reference frequency with Ky assumed infinite and the other of finding 
the significant rotor speeds referred to the other support frequency 
| /%/M as reference frequency with K x assumed zero. 

5. The analysis of a four-blade counterrotating rotor system in 
which the rotors cross along the principal stiffness axes of the rotor 
supports leads to the same equations as those considered for the special 
case of Ky = and the 'stability properties are given by the investi- 
gation of that special case. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., August 24, 1956. 
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APPENDIX A 
SYMBOLS 


a 


radial position of vertical hinge 


a^Q, so f° r th 

A l> B l> c l 
®Z > * 

b 

^* n< ^ 

. so forth 
D 

F(c%) = lim F n (c%) 
n >00 

Fn ( ^) = ~ 

j, n, s 
J(afe) 

k = 2A 5 - 1 
K(oe./co) 

K(0), K(l), K(l/2) 

Kl, K 2 


elements of determinant defined in equa- 
tion (l6) 

Fourier coefficients in equations (12) 

complex conjugates of A i, B and C 
respectively 

distance from vertical hinge to center of 
mass of blade 

minors of the determinate ^((%) 

time -derivative operator, — 


integers 

function of F(ofe) defined by equation ( 27 ) 

function of ofe/o) defined in equation ( 29 ) 

functions defined by equations ( 30 ) 
functions defined by equations (35) 


K x , Ky 


spring constants of the rotor supports along 
the X- and Y-directions, respectively 
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average stiffness of rotor supports. 


Ky + Kx 


K v - K* 


spring constant of blade self-centering spring 


effective mass of rotor supports 


mass of rotor blade 


total mass of two-blade rotor system, m + 


P(t), Q(t), R(t) periodic functions defined in equations (ll) 


P(t), Q(t), R(t) 


complex conjugates of P(t) , Q(t), and R(t), 

respectively 

constant defined in equations (33) 

radius of gyration of blade about center of mass 


x r> Yr 


kinetic energy 

potential energy 

deflection of rotor hub measured in fixed 
X,Y-coordinate system 

deflection of rotor hub measured in rotating 
X,Y- coordinate system 

fixed rectangular coordinate axes taken parallel to 
the principal stiffness directions of the rotor 
hub 

complex position coordinate of the rotor hub in 

rectangular coordinate system rotating with angular 
velocity, x r + iy r 


complex conjugate of z, x r - iyr 


jj 
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Pi, P2 angular hinge deflections of rotor blades, respectively 

A discriminant of cubic equation x 3 + bx 2 + cx + d = 0 

A(c%) determinant of infinite order defined by equation (16) 

£*i(cjc&) determinant of order 6n - 3 formed from A(c%) 

e o = ^(Pi + P2) 

01 = |(Pl - P2) 


a 




^i> 'ii> ^0 


blade variables for counterrotating rotor 



, . ^b ’ 
mass ratio, 

M 


0 ) 

<% 


constant angular velocity of rotor 

characteristic exponent or natural frequency of rotor 
system as viewed in coordinates rotating at angular 
velocity u) 

principal values of <% 
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APPENDIX B 

THE GENERAL EQUATIONS OF MOTION FOR TWO-BLADE ROTORS 
By George W. Brooks 


General Case 

The present report, although it generalizes the theory of ground 
vibrations of two-blade rotors by treating the case of anisotropic 
supports, does not consider the effects of damping. Reference 2 treated 
the case of isotropic supports or equal pylon stiffness in all horizontal 
directions and included damping. The more general case is one which 
involves the treatment of the two-blade rotor with anisotropic supports 
and damping. The equations of motion for this case, derived in the com- 
plex coordinate notation of the present report, are given here. 

The equations for the kinetic and potential energies of the rotor 
system are given by equations (5) and (6) of the text as follows: 


= ^(z + ioiz)(z - icuz) + 
2 


m b 


(z - icoz)(i0i - cu 9 i) - (z + iooz)(i§i + 061) + 


( x + + §l2 



(Bl) 


v=|z, + !i(e 0 2 + ei 2) 


£K( z 2 e 2icot + z 2 e -2ioJt 

4 ' 


(B 2 ) 


Two types of damping of the rotor system are assumed to exist: 

(l) damping in the rotor supports which is proportional to the transla- 
tional velocity of the rotor hub and (2) damping of the rotation of the 
blades about the drag hinges. The damping dissipation function is then 


F = — 


2 L 


zz + a2zz + iu>(zz - zz) 


B, 




+ 


(B3) 
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where B is the damping force per unit velocity of rotor-hub displace- 
ments and Bp is the damping torque per unit angular velocity of blade 
motion about the drag hinge. 

By using the Lagrangian form of the equations of motion 


dT + dF + dV = 0 (B4) 

dt \dx ) dx dx dx 

and usiug as variables 8q, z, z, and 0^, respectively, the following 
equations of motion for the rotor system axe obtained: 


(d 2 + DAp 


cd 2 A]_ 


+ A2)0 o - 0 


(B5) 


(D + ioo) 2 + A(D + io>) 



— ze“2iojt + i^(p + iao) 2 0i = 0 

M 


(B6) 


— ze 2ia5fc + 
M 


(D - ia>) 2 + A(D - ioo) + 


z - in(D - itu) 0 X = 0 (B7) 


(D + ia)) 2 z - (D - icu) 2 z + 2i(l + 


^j^D 2 + DAp + A^cd 2 + A 2 J 0i = 0 (B8) 


where two additional combinations of parameters have been introduced, 
namely. 



A = 


B 

M 


and 
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Special Case - Zero Damping 

If damping is neglected, Ap = A = 0 and the equations of motion, 
(B5) to (b 8), reduce to 


(D 2 + o?Ai + A2)9 q = 0 


(B9) 


(D + i<u) 2 + ^ 
M 


- ^ ze~ 2ia}t + ip.(D + ico) 2 ©! = 0 (BIO) 


_ ^ ze 2ia>t + 
M 


(D - io)) 2 + - 
M 


z - in(D - io)) 2 e 1 = 0 (Bll) 


(D + icn) 2 z - (D - io>) 2 z + 2ifl + — HD 2 + A]o£ + A2)0l = 0 (B12) 


which are the equations of motion for the case treated in the present 
report, that is, equations (7 ) 'to ( 10 ). 


Special Case of Isotropic Supports 
If K x = Ky, then 2K = 0 and equations (B5) to (B8) reduce to 


(D 2 + DAp + + A2)0Q = 0 


(B13) 


(D + ioo) 2 + A(D + icn) + — 

M 


K 


z + ip(D + ico) 2 0i = 0 (B14) 


(D - io>) 2 + A(D - io>) + & 

M 


z - ip(D - i(u) 2 0! = 0 (B15) 


(D + icn)^z - (D - io>) 2 z + 2i (l + j(D 2 + DAp + A]_o^ + A2)0q = 0 (Bl6) 

h 2 / 


r2\ 
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Equation (B13) involves only the in-phase motion of the blades, can 
be solved separately, and will not be further considered here. If z 
and z are expressed in the equivalent notations x + iy and x - iy, 
respectively, and equations (Bl4) and (B15) are added together, equa- 
tion (l) of reference 2 is obtained. Equation (2) of reference 2 is 
obtained by direct substitution and equation (3) is obtained by subtrac- 
tion of equation (B15) from equation (Bl4) . Thus the equations of motion 
for the case treated in reference 2 as well as for the case treated in 
the present report are shown to be special cases of the general equations 
of motion, equations (B5) to (B8), for a rotor mounted on anisotropic 
supports with hub and blade -hinge damping. 
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' APPENDIX C 

DERIVATION OF THE CONVERGENCE FACTOR FOR K(0) 

A convergence factor for K(0) is found by finding a simpler func- 
tion Gq that changes with n in nearly the same way as K(0) n . Then, 

if G denotes lim the expression 


K(0) n G 

G n 


for a given value of n is a better approximation to K(0) than is 
K(0) n alone. A suitable form for Gn is found from a study of the 
behavior of K(o) n for large values of n. 

The behavior of K(o) n is studied by first observing the behavior 
of C n (0) and Ln(0) for large values of n and then inferring the 
behavior of K(0) n from the appropriate relation of equations (52). In 
the discussion of equations (22) it was shown that, as n becomes infi- 
nite, the ratio C n+ i^C n approaches the value k (eq. (22a)). A closer 

approximation to the value of this ratio can be written as 


= k 4 P _Q_\ 

C n (0) \ n kn 2j 


where P and Q are constants to be determined. 
Equations (19) become, for c% = 0, 


C n +l(0) = 

Dn+l( 0 ) = 

E n+l(°) = 


-1 + 


k/m 


-1 + 


(2n + 1 ) d a£ 
Aju£ + Aq 


Dn(0) - A3 


f 2n + 

[ 2n i 


E n( 0 ) 


-1 + 


(2n) 2 o£ 

k/m 


(2n - 1) a? 


En(0) - A 3 

c n (o) 


f 2n - 1 
l 2n . 


Cn(0) 


(Cl) 


(C2) 
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where the second term in the equation for En(a&) in equations (19) has 

been neglected as being of higher order in powers of l/n than the terms 
retained. Eliminating E(0) and D(0) from the equations (C2) results 
in 


Cn+l(C) _ 

1 1 K/k 


An I (JC? + Ao 

1 1 

1 , K M... 

c n (o) 

(2n + l) 2 o£ 


(2n) 2 o£ 

(2n - l) 2 o£ 


(2n - if 

\ , K / M 

+ fa + A 2 

\ 

fit K /M 

V ) 

\ 

(2n + 1) 2 o£ 

\ 2n ) 

(2n - l) 2 o£ 

/ 


Upon expanding this expression into powers of l/n and retaining terms 
up to and including those in l/n^ 


c n+ i(o) 

c n (0) 


-1 + 2A 3 - 

' (2n) 


2 | + Aio^ + A2 




Comparing equation (C3) with equation (Cl) ■ results in 


P = 0 


(C3) 


2 &(l - A 5 ) + Aqu2 + A 2 + 2Aju2 
o^(l - 2Aj) 

Similarly 


Ln+l(0) 

4n(0) 


Q 

(2n) 2 


Therefore, from equations (32) 


K(0) n+ i 



K(0) n 
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The right-hand side of equation (C4 ) , which is seen to be of the form 
G/G n , is a convergence factor. This convergence factor is the one used 

in equations (33)- 



NACA TN 3844 


113 


APPENDIX D 


MATHEMATICAL ANALYSIS FOR THE CASE OF Ky = ~ 

In terms of the variables 0 q> 9^, and x, the expressions for the 

kinetic and potential energy are 


T = ^ x 2 + m b 


rW +h2 j 


- o£ ^(e 0 2 + ©I 2 ] 

b 


2x(§]_ sin cut + o39i cos cut) 


= — X 2 + ^ (q 0 2 + Qi 2 ] 

2 b 2\ / 


The three equations of motion are 


0q + 9q = 0 


.. K x 

X + X - LI 

M 


- — (©1 sin cut) = 0 


dt 


2 ' 


(Dl) 


(D2) 


- — - x sin cut + 0q + ^cjcRAi + A2^0i = 0 


1 + *§ 

b 2 


(D3) 


Equation (Dl) is identical with equation (7)- Equations (D2) and (D3) 
constitute a system of two linear second-order differential equations 
with periodic coefficients. 
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ll4 


Equations (D2) and (D3) axe satisfied by solutions of the form 


X = 2Z A l e ± ^ +la ^ t 
l=-CO 


— i fo^a+( ^+l)o)j t 
01 = Z B i+i e 

l =-00 


(d4) 


where l takes on all odd integral values and Ai and are con- 

stants to be determined along with the two principal values of 

and ^2)' cons 4ants A^ and B , in equations (d 4) are, 

of course, different from those in equations (12). 

Combining equations (D4) with equations (D2) and (D3) and setting 
the coefficients of the various exponential time factors equal to zero 
give 


-1 + 


K x M 


2i 1 + — 

V b2 y 


(<^a + Ic 0) 

c% + (l - 2)cu 


A l - ^("B l-l + B l+1 ) = 0 


<% + (l - 1) 


(JO 


h-2 + 


C% + l(JD 


<% + (l - 1) 


A Z S' + 




-1 + 


a?A]_ + A2 


(c% + (l - 1)0^' 


> B l-1 = 0 
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The determinant al equation is then equal to 



• • • a - 3,-3 a - 3,-2 0 

• ' * a -2,-3 a -2,-2 a -2,-l 

• • • . 0 a -l, -2 a -l,-l 

. . . 0 0 ao^_i 

• • • 0 0 0 

• • • 0 0 0 

• • • 0 0 0 


0 

0 

0 

0 

0 

0 

0 

0 

a -l,0 

0 

0 

0 

®0,0 

80,1 

0 

0 

a l,0 

a l,l 

a l,2 

0 

0 

3-2,1 

a 2,2 

a 2,3 

0 

0 

a 3,2 

a 3,3 


= A(a^) = 0 


where 


-3,-3 


.-i + 


(c% - 3a))‘ 


A-icu? + Ap 
30,0 = -1 + ^- 72 — 




a 3 

a -5,-2 = — 


a 0,l = 


f <% + (o' 

V ^ J 


- 3(0 ^ 

i " 2, “ 5 Vafe - 2o)J 


_ a 3 
a i,o - — 


< 
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- 2,-2 = -1 + 


(c% - 2(d) ‘ 


K x M 


a l 1 = "I + ^ 

((% + Cd) 


l -2,-l = 


<% - 2 (JD 


- A 3 

a -l,-2 = ^ 


/(% + CD 

3,2,1 \c% + 2(d 


K x /m 


‘-1,-1 - - x o 

(<% - cd) 


At + Ao 

a 2 2 =- 1 + p 

'(<% + 2a>) 


a 3 

a -l,0 = — 


<% + 3a> 

(<% + 2co) ‘ 


^j-l 


C% - CD 


'3,2 - — 


a 3,3 “ " 1 + 


Kx/M 

t 

(c% + 3 cd)' 


o£A]_ + Ap 

The term -1 + will be taken as the origin of’ the determinant. 

a* 2 

Define ^(a^.) as the determinant of order 4n-l formed by taking 
a square array from A(oca) centered on the origin. Then 


A(c%) = lim ^a(c%) 
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Define auxiliary determinants from ^(a&) as follows: 

Cn(oa) determinant ' of order 2n-l consisting of the terms to the 
right of and below the origin term 


Dn(ofe) 

determinant of order 2n-2 
its last row and column 

obtained 

from 

Cn(<%) 

by 

omitting 

Mn(ofe) 

determinant 
its first 

of order 2n-2 
row and column 

obtained 

from 

Cn(<m) 

t»y 

omitting 

Nn((«a) 

determinant 
its first 

of order 2n-3 
row and column 

obtained 

from 

Dn(ufe) 

by 

omitting 


The determinants C n (a^) and Dn(c%) satisfy the following recur- 
rence relations: 


l “1 + 


Kx/M 


[<% + (2n - 1 )cd] £ 


l Dn(ofe) - 


A 5 

2 


DnC^h) = S 


+ (2n - 1) 


-i2 


0) 


<% + (2n - 2 )od 

+ A 2 


Cn-l(<%.) 


-1 + 


jt% + (2n - 2) os 


> Cn-l(oJa) ~ 


^3 

2 


03 a + (2n - 3) 


CD 


o^ + (2n - 2) 


CD 


Dn-lC^a) 


(D5) 


The recurrence relations (eqs. (D5)) are also satisfied by • M n (<%) and 

Nn(afe)? with M and N replacing C and D, respectively. The initial 
values are 

kJm 


Cl(tha) = -1 + 


(c% + 03 ) ‘ 
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= 


^2(0^) = -1 + 


K X J 

[m 

cjg^Ai + A2 

_ % 

2 

+ <•> \ 

X - 

(<% ■ 

v cd ) 2 

(c% + 2o>) 2 
— J 

2 ' 

y% + 2 &) 


u^A^ + Ag 
(<% + 2 cjd) £ 


M2(cwa) = 


+ A2 

. K x) 

M 

_ ^3 

+ 3o)\ 

2 

(<% + 2<n) 

(<% + 

3m) 2 

2 1 

Uda + 203/ 


Expanding ^n(ofe) in terms of the elements of the column containing 
the origin gives 


An(aa) = -1 + 


oc?A-j_ + aX 


At: 


<% 2 j 

\2 


Cn(-cta)C n (c%) - 


— -) C n (o^)M n (-t%) - ( — + — ) C n (-a^)Mn(afa) 




o& 


As n becomes infinite, or as <% becomes infinite, the recurrence 
relations (eqs. (D5)) approach 


Ax 

Cn = -Dn - ~ Cn-1 


Ax 

Dn = -Cn-1 ~ “ ^n-1 


(D6) 


Equations (d 6) are satisfied by a solution of the form 

C n = Ck n 


Dn = -c(k + -|jk n_1 


(DT) 
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t 


vhere k satisfies the equation 


A* 

- (1 - A3)k + -j- = 


0 


(D8) 


The larger root of equation (D8) will be denoted by k and the smaller 
root by kp. Although the complete solution of equation (d 6) is of the 
form 


C n = Ck n + Ck^ 

for large values of n the term in kq becomes negligible compared with 
the term in k. 

With the same values for k, M n and N n will have solutions similar 
to those of equations (D7). Thus as n becomes infinite will 

vary as the quantity k 2n . 

Define the function 


F(aca) = lim Fn(ofe) 
n — >oo 



Ai(<ha) 

k 2n 


The function F(oc&) is periodic in (% of period 2oo, has roots 
±(o)ai - 2so)^ i(oca2 - 2sco), and second-order poles at . (oca ± sod) for 
all integral values of s. Furthermore, . F(ooa) approaches the limit 


E = lim F(oo£l) 
c 


lim F n (oo) 

n — ^oo 


. c 2h(°°) k 2 

= lim = 

n — k.2n k] 2 - k 2 

as becomes infinite in a direction other than along the real axis. 
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Form the function 



(D9) 


The function J(cca) is an analytic function of everywhere. Hence, 

by Liouville's theorem, J(u&) is a constant. By letting — >+ ® 

along the imaginary axis, it is seen that J(oja) = -4E. 


Substituting -4E for j(oa) into equation (D9) results in 


11m = F(<%) 

n — k 211 



Introducing K functions defined similarly to those' used for the 
case Ky = 0 gives 


. . 2 /”V 

- a — sinM 


Ki = sin^l- 


\ i s m t 

\ 2a) 1 \ 2a) 


lim 


s “K^) f(c%) 


-4e 


K2 - cos 


2 |^i 

V 2m 




cos £ 


\ 


lim 
<% — 


sin c l - — 
V2co > 


2a> 
|F(afe) 


-4e 




(DIO) 
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Carrying out the limit processes indicated, in equations (DIO) gives 


Kq = lim K in 
n — >oo 


= lim 
n — >oo 


+ A2)c n 2 (0) - A3* 2 o£c n (0)N n (0) 


-4a^C2n(“) 


K2 = lim K2- 


n — >oo 


lim 
n >00 


■n 


* 2 ^ 


-4a> 2 C2n(«>) 


Finally, upon introducing appropriate, convergence factors, the 
following quantities needed in equations (36) are obtained: 


Kq_ = lim 

n — >00 


K X sta 2 ('’f' 
ln \2 ) 


rt 2 R 2 iV /, R 2 


TT(i- 

0=1 \ 4j 


*2 


K 2 = lim ^ 
n — >00 


K 2n “ s2 (f ) 


n 


IT 

0=1 


1 - 


R 2 

(2J - 1) : 


where 




-jjj* + u£Aj_ + Ag + oc^^l + k]_ - kj 


o£(k - ki) 
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x Undeflected 
hub -position 


Figure 1. - Simplified mechanical system representing rotor 



Figure 2.- Chart for obtaining principal value of from yalues of 

Ki and K 2 - K y = 0 or K v = 


00 . 
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Figure 8.- Response frequencies of a two-blade rotor on symmetric supports 
for Aj = 0.1, A2 = 0, Aj = 0.1, and Ky = K^.. 
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Figure 9.- Chart giving position of the main instability range for K, 
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